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SUMMARY  PAGE 


THE  PROBLEM 

Digitized  photo  data  acquired  during  Impact  experiments  must  be  scaled  and 
converted  to  a  format  suitable  for  analysis.  Once  converted,  these  data 
require  comprehensive  graphical  presentation  for  efficient  Interpretation. 
This  report  presents  a  detailed  description  of  the  software  developed  to 
accomplish  such  tasks  In  both  a  production  and  an  Interactive  environment. 

FINDINGS 

Procedures  utilizing  the  design  presented  have  been  Instituted  at  NBDL  and 
found  to  be  effective. 

RECOMMENDATIONS 

In  an  environment  requiring  data  reduction  processing  of  photo  data,  the 
design  presented  herein  should  be  considered. 
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A  METHOD  OF  PHOTO  DATA  REDUCTION,  WITH  DESIGN  CONSIDERATIONS 
FOR  THE  NOVA  800*  AND  UNIVAC  1100/83®  COMPUTERS 


INTRODUCTION 

The  Naval  Biodynamics  Laboratory  (NBDL),  located  In  New  Orleans,  LA,  is  an 
internationally  recognized  laboratory  which  performs  experimental  research  to 
determine  the  effects  of  aircraft  crashes,  ship  motion,  vibration,  aircraft 
ejection  and  parachute  opening  forces  on  the  health  and  performance  of  Navy 
personnel.  On-going  research  programs  use  high  speed  Instrumentation  cameras 
to  record  the  motion  of  test  subjects  during  biodynamic  experiments.  The 
films  are  digitized  and  the  3-dlmenslonal  motion  is  reconstructed  and 
analyzed. 

The  procedures  and  programs  used  to  convert  the  digitized  photo  data  Into  a 
format  compatible  with  a  large-scale  general  purpose  computer  (UNIVAC  1100®), 
to  plot  the  data  for  identification  of  data  errors,  and  to  prepare  the  data 
for  subsequent  processing  have  recently  undergone  major  revision.  The 
objectives  of  the  revision  were  to: 

1.  Modify  existing  software  to  be  compatible  with  announced  changes  In 
the  UNIVAC  operating  system. 

2.  Use  new  graphics  support  capabilities  to  Improve  error  detection 
procedures,  reduce  turn-around  time  of  analysis,  and  produce  output  compatible 
with  all  available  graphics  devices. 

3.  Improve  efficiency  of  operations  personnel. 

4.  Standardize  data  formats  and  I/O  access  methods  to  be  compatible  with 
existing  archival  data. 

5.  Develop  a  modular  organization  of  software  which  provides  greater 
operational  flexibility. 

This  report  documents  the  new  procedures  and  programs. 

DESCRIPTION  OF  SOFTWARE 

Preparation  of  digitized  photo  data  for  subsequent  reduction  and  analysis 
requires  a  process  consisting  of  several  tasks  performed  under  computer  con¬ 
trol.  Software  supporting  each  task  conforms  to  the  principles  of  modular 
design  so  that  task  sequences  may  be  selected  to  meet  specific  needs  of  the 
data  being  processed. 

In  order  to  process  digitized  photo  data  through  a  data  reduction  flow,  major 
tasks  are  selected  and  merged  In  a  single  runstream.  The  three  available  pro¬ 
grams  are: 

1.  Convert  photo  data  from  PDS  (Photo  Digitizing  System)  NOVA  800® 
format  to  UNIVAC  1100/83®  format  for  subsequent  analysis. 


» 


2.  Produce  film  frame  plots  and  x-y  trajectory  plots  of  target  data  in 
order  to  provide  graphical  presentation  for  error  analysis. 

3.  Convert  data  digitized  manually  into  the  format  of  data  tracked  under 
computer  control . 

Each  of  these  components  is  documented  in  detail  in  the  body  of  this  presenta¬ 
tion.  It  is  assumed  that  the  reader  is  familiar  with  the  operating  system  of 
the  UNIVAC  1100/83*  computer  (see  reference  3)  and  with  the  UNIYAC  ASCII 
FORTRAN  programming  language  (see  reference  4).  Although  the  software  pre¬ 
sented  herein  is  intended  for  specific  computers,  plotting  devices  and  soft¬ 
ware  support  packages,  it  has  been  designed  to  minimize  the  conversion  effort 
for  other  systems. 

1.  Photo  Data  Conversion  Program 

A.  Introduction.  This  section  describes  a  UNIVAC  1100/83®  program  which 
converts  photo  data  from  PDS  (Photo  Digitizing  System)  data  tapes  to  UNIVAC 
1100/83®  format.  In  converting  data  to  a  format  compatible  with  a  large-scale 
general  purpose  computer  (UNIYAC  1100/83®),  the  photo  data  conversion  program 
provides  a  suitable  environment  for  efficient  data  reduction  and  graphical 
presentation. 

B.  Discussion.  Input  to  the  conversion  program  is  a  PDS  data  tape.  It 
contains  digitized  photographic  data  from  at  least  one  experiment,  or  run. 

For  a  particular  run,  as  many  as  three  cameras  are  used  to  photograph  targets 
on  anatomical  mounts  strapped  to  the  experimental  subject.  Each  reel  of  film 
generated  from  these  cameras  during  the  run  is  digitized  (on  the  NOVA  800® 
system),  with  the  output  being  written  to  a  PDS  data  tape.  Data  for  several 
runs  may  be  written  to  a  single  PDS  tape.  This  tape  contains  four  record 
types.  The  first  is  a  header  record,  which  contains  identification  informa¬ 
tion  for  film  from  a  particular  camera  used  during  a  run.  This  record 
includes  such  Items  as  run  number,  PDS  record  type,  camera  site  number  (1,2 
or  3),  jullan  date  the  film  was  digitized  and  time-of-day  for  the  first  film 
frame  (all  character  data  on  the  PDS  tape  is  in  ASCII  format  and  numeric  data 
is  in  twos  complement  integer  format).  The  second  record  type  is  the  time 
record,  which  identifies  time-of-day  for  each  film  frame,  or  data  point.  The 
time  data  values  originate  from  LED  displays  on  the  film  frames  themselves. 
There  are  250  data  points  on  a  time  record  (the  first  point  being  the  time  at 
first  motion  of  the  subject).  Identification  information  for  the  reel  of  film 
the  time  record  describes  is  also  Included  in  the  record,  i.e.,  run  number, 

PDS  record  type  and  camera  site  number.  The  final  two  record  types  are  the  X 
and  Y  data  records,  respectively.  Each  contains  250  data  points  measuring 
film  plane  displacement  (Inches)  from  a  stationary  origin  on  the  film.  An  X 
and  Y  data  record  exists  for  every  photo  target  tracked  during  the  run.  Each 
data  record  contains  (in  addition  to  displacement  values)  run  number,  PDS 
record  type,  camera  site  number  and  a  target  number.  The  target  number  con¬ 
tains  Information  describing  the  anatomical  mount  type  on  which  It  is  located. 
For  Instance,  there  are  head  mounts,  neck  mounts,  mouth  mounts  and  pelvic 
mounts  -  all  designed  to  track  motion  of  a  particular  segment  of  the  body. 


Data  on  a  PDS  output  tape  Is  organized  by  reels  of  film.  That  Is,  all  data 
for  a  particular  reel  of  film  (with  Its  unique  header,  time  and  data  records) 
are  grouped  together  Into  a  file.  The  order  In  which  the  records  appear  Is 
always  the  same:  header  record,  time  record,  target  1  (x  values),  target  1  (y 
values)  •••  target  N  (x  values),  target  N  (y  values),  where  N  =  number  targets 
digitized  on  film.  In  some  experiments  more  than  the  standard  250  frames  are 
digitized.  In  this  Instance  the  data  records  for  the  extra  frames  (always  a 
multiple  of  250)  are  written  to  the  PDS  tape  Immediately  after  the  final  data 
record  containing  data  for  the  previous  250  frames.  An  end-of-file  separates 
each  group  of  data  records  for  a  particular  set  of  250  frames.  A  double  end- 
of-file  signals  the  end  of  data  on  the  tape. 

The  conversion  program  reads  the  NOVA  800®  input  tape  and  stores  It  in  an 
Input  buffer.  The  PDS  characters  are  In  8-bit  ASCII  format  and  the  numeric 
data  Is  either  In  16-bit  or  32-bit  integer  format.  The  conversion  program 
determines  the  format  of  the  Information  to  be  converted  from  the  format 
specifications  of  the  PDS  Input  tape  (see  Figures  lb-ld).  All  character  data 
Is  converted  from  8-bit  NOVA  ASCII  to  9-bit  UNIVAC  ASCII  format.  Numeric  data 
is  either  converted  to  36-bit  UNIVAC  floating  point  format  or  36-bit  UNIVAC 
Integer  format,  depending  on  the  format  type  of  the  Information  being  con¬ 
verted.  For  Instance,  PDS  record  type  Is  converted  from  NOVA  16-bit  Integer 
to  UNIVAC  36-bit  Integer,  while  displacement  (data)  values  are  converted  from 
NOVA  32-bit  Integer  to  UNIVAC  36-bit  floating  point  format.  All  converted 
displacement  values  are  then  scaled  from  inches  to  meters.  This  Is  done  to 
provide  data  compatibility  with  a  group  of  data  analysis  programs  collectively 
called  EASYFL0W*. 

The  converted  and  scaled  photo  data  must  be  output  In  a  standard  file  format 
suitable  to  subsequent  processing;  In  addition  this  data  must  be  output  In 
FIELDATA  format,  which  Is  the  NBDL  archival  standard.  This  second  requirement 
presents  difficulties  in  an  ASCII  FORTRAN  programming  environment.  In  order 
to  provide  for  both  needs  an  1/0  access  package  (see  Reference  1)  has  been 
developed  to  (1)  provide  standard  1/0  access  methods  for  an  established  photo 
data  file  format  and  (2)  maintain  data  format  compatibility  with  the  NBDL 
archives.  The  1/0  package  provides  directory  maintenance  for  the  user.  The 
conversion  program  utilizes  this  capability  to  output  the  standardized  direct 
access  record.  This  directory  contains  unique  run  Identification  numbers  (for 
all  runs  on  the  output  tape)  and  accompanying  sector  locations  for  both  header 
and  time  records  for  the  run.  The  directory  provides  means  for  rapid  direct 
access  of  photo  data  for  a  given  run  number. 


*  EASYFL0W  Is  an  NBDL  term  representing  a  data  reduction  and  analysis  system 
requiring  Input  data  of  a  specific  format.  EASYFLOW  uses  converted  photo 
data  to  derive  position  coordinates  In  3-dlmenslonal  space  of  targets 
affixed  to  anatomical  mounts  (which  are  strapped  on  to  the  experimental 
subject).  These  data  are  subsequently  used  to  compute  and  plot  displace¬ 
ment,  linear  (and  angular)  velocity  and  linear  (and  angular)  acceleration 
of  the  subjects. 


There  are  four  record  types,  in  addition  to  the  directory,  which  the  conver¬ 
sion  program  outputs  to  a  UNIVAC  mass  storage  file.  The  first  of  these  is  the 
header  record.  There  exists  a  unique  header  record  for  each  camera  site  used 
in  a  run.  There  are  no  more  than  three  and  no  less  than  two  cameras  used  In  a 
run,  each  assigned  a  site  Identification  number  of  1,  2  or  3.  This  site  ID 
and  the  run  number  In  concert  make  the  header  record  unique.  A  keyword  is 
stored  In  the  header  record  along  with  Identification  information  such  as  the 
run  number  and  site  ID.  The  keyword  identifies  the  record  as  a  header.  The 
second  type  of  output  record  is  the  time  record.  It  contains  the  time-of-day 
information  for  each  of  250  film  frames.  In  addition  the  time  record  includes 
the  run  number,  the  site  ID  and  a  keyword  for  identification.  The  third 
record  type  is  a  photo  data  record.  It  contains  displacement  values  (from  a 
fixed  point  on  the  film  plane)  in  meters  of  a  particular  photo  target  on  the 
reel  of  film.  Also  included  are  the  run  number  and  site  ID  of  the  film,  along 
with  a  keyword  identifying  the  record  type.  This  keyword  also  contains  infor¬ 
mation  identifying  the  anatomical  mount  (on  which  the  photo  target  resides) 
and  the  axis  type  of  the  data  (X  or  Y).  The  combination  of  header  record, 
time  record  and  data  records  (for  all  targets  on  the  reel  of  film)  is  written 
to  disk,  in  the  order  listed,  by  the  conversion  program.  Data  for  each  reel 
of  film  with  a  camera  site  ID  of  1  are  grouped  together  on  the  output  file  and 
separated  by  double  ends-of-file*.  An  additional  end-of-flle  is  written  after 
the  final  data  record  of  the  final  reel  of  film  for  site  ID  1.  This  signifies 
end-of-data.  However,  the  direct  access  table  follows  this  final  end-of-flle. 
It  is  located  at  the  end  of  the  file  so  as  not  to  limit  its  size.  In  order  to 
provide  access  to  the  directory  a  "locator  table"  is  written  as  the  initial 
output  record.  It  contains  the  sector  address  of  the  directory. 

The  above  output  file  structure  is  duplicated  for  camera  site  2  (and  3  if 
necessary).  A  separate  UNIVAC  disk  file  is  thus  created  for  each  site  ID. 

Each  of  these  files  is  in  a  standardized  format,  and  serves  as  input  to  a 
UNIVAC  runstream  designed  to  write  them  to  a  UNIVAC  1100/83®  data  tape. 

The  conversion  program  is  designed  to  correct  three  common  situations  in  which 
PDS  operator  error  has  occurred.  These  three  special  cases  are: 

(1)  The  run  number  of  a  header  record  for  a  reel  of  film  is  entered 

incorrectly  by  the  PDS  operator  and  written  to  the  PDS  data  tape. 

(2)  A  reel  of  film  is  re-dlgltlzed  by  the  PDS  operator  and  output  to  the 

end  of  the  same  PDS  data  tape  as  data  from  the  original  reel  of  film. 

(3)  The  camera  site  ID  of  a  header  record  for  a  reel  of  film  is  entered 

incorrectly  by  the  PDS  operator  and  written  to  the  PDS  data  tape. 


*  If  more  than  250  frames  of  data  were  digitized  on  a  reel  of  film,  each 
group  of  data  records  for  the  extra  frames  (always  a  multiple  of  250)  is 
written  in  succession  and  separated  by  single  end-of-flle. 


In  cases  (1)  and  (3)  the  conversion  program  uses  correction  fields  In  card 
Image  Input  In  order  to  scan  the  header  record  for  the  Incorrect  entry  and 
then  correct  It.  This  Is  done  prior  to  writing  the  header  to  disk.  In  case 

(2)  the  conversion  program  scans  for  the  first  occurrence  of  the  data  file  for 
a  reel  of  film.  Once  found,  this  data  file  (header,  time  and  data  records)  Is 
bypassed  and  thus  not  converted  and  output.  However,  the  second  occurrence  of 
a  data  file  for  the  reel  of  film  is  converted  and  output  in  the  usual  manner. 
The  photo  conversion  program  accomplishes  the  following  objectives: 

(1)  Converts  photo  data  to  a  format  compatible  with  a  large-scale  general 
purpose  computer  (UNIVAC  1100/83*) . 

(2)  Scales  displacement  data  values  from  Inches  to  meters  In  order  to 
provide  data  compatibility  for  subsequent  processing. 

(3)  Provides  for  correction  of  common  PDS  operator  errors. 

(4)  Accomplishes  output  of  data  with  standardized  I/O  access  methods. 

(5)  Maintains  data  format  compatibility  with  the  NBDL  archives. 

C.  Main  Program  and  Subroutine  Descriptions: 

(1)  MAIN  (see  Appendix  lb  for  listing)  -  This  is  the  driver  program 
which  controls  all  program  operations.  All  correction  cards  are  read  by  this 
program  and  used  to  rectify  the  three  special  cases  earlier  noted.  MAIN  then 
calls  subroutine  CONVT  In  order  to  read  and  convert  the  first  PDS  Input  record 
(see  Figures  le-lg  for  Input  record  formats).  The  rest  of  the  input  records 
(non-header)  for  the  run  are  read  and  converted  by  virtue  of  calls  to  sub¬ 
routine  CNVT1.  Here  the  converted  records  are  written  to  mass  storage.  Once 
a  run  has  finished  processing,  control  returns  to  MAIN  and  the  process  repeats 
Itself  until  all  runs  have  been  processed.  At  this  point  a  triple  end- 
of-flle,  denoting  end-of-data.  Is  written  to  the  output  file. 

(2)  CNVT1  (see  Appendix  lc  for  listing)  -  This  subroutine  controls  the 
conversion  and  output  of  all  data  for  each  camera.  The  main  program  reads  the 
header  records  and  transfers  control  to  this  subroutine  for  reading  and  con¬ 
version  of  the  position  and  time  records.  The  position  data  Is  converted  from 
Inches  to  meters*.  Each  header,  time  and  position  record  Is  written  to  mass 
storage  by  a  call  to  the  photo  10  subroutine  PUTRCD  (see  Appendix  4c).  Data 
for  each  site  ID  Is  written  to  a  separate  mass  storage  unit  (see  Figure  2a). 

As  mentioned  earlier,  the  necessary  data  directories  are  maintained  and 
written  by  the  photo  10  subroutines.  After  processing  an  entire  run,  CNVT1 
writes  a  double  end-of-flle  for  each  output  unit  and  returns  control  to  the 
calling  program.** 


*  During  the  scaling  process  all  data  values  of  0.0  (points  not  digitized) 
are  set  to  999.0 

**  If  more  than  250  frames  were  digitized  for  each  target  In  a  run,  only  a 
single  EOF  (End-of-flle)  Is  written  after  the  final  data  record  for  each 
set  of  250  frames. 


(3)  CONYT  (see  Appendix  Id  for  listing)  -  This  subroutine  is  called  by 
CNVT1  to  perform  the  actual  conversion  of  data  from  PDS  format  to  UNIVAC  1100/ 
83*  format.  It  is  called  after  each  read  of  a  PDS  input  record.  All 
character  data  is  converted  from  8-bit  ASCII  to  9-bit  ASCII  format.  Numeric 
data  is  converted  from  32-bit  binary  integer  to  UNIVAC  1100/83®  36-bit  binary 
floating  point  format  (see  Figures  2b-2f  for  output  record  formats). 

In  this  subroutine  the  sign  bit  is  extended  for  negative  numbers;  these  nega¬ 
tive  values  are  data  entered  manually  by  the  PDS  operator  by  using  the  X,Y 
crosshair  input. 

(4)  SPHAVG  (see  Appendix  le  for  listing)  -  This  subroutine  calculates 
the  X  and  Y  position  of  a  stationary  reference  target  relative  to  the  film 
sprocket  hole  for  each  camera  used  in  a  run.  Sprocket  hole  values  for  five 
consecutive  frames  are  averaged,  the  result  being  subtracted  from  the  averaged 
sled  (chair)  displacements  for  the  same  five  frames  in  order  to  produce  the 
final  value.  The  X  and  Y  sprocket  hole  averages  are  printed  for  each  camera 
and  in  those  cases  where  there  is  no  sprocket  hole  information  on  the  PDS 
tape,  the  values  can  be  read  from  the  film  by  the  PDS  operator  and  calculated 
by  hand.  Execution  of  EASYFLOW  is  expedited  by  providing  sprocket  hole 
averages  at  the  time  of  conversion.  This  information  is  required  for  an 
internal  consistency  check  on  the  data  collection  and  system. 

(5)  ASC8T9  (see  Appendix  If  for  listing)  -  This  subroutine  converts  a 
string  of  8-bit  ASCII  characters  to  a  UNIVAC  1100/83®  compatible  9-bit  ASCII 
character  string.  All  character  data  read  from  the  PDS  input  tape  is  con¬ 
verted  by  the  subroutine.  The  high  order  bit  of  a  9-bit  ASCII  character  is 
always  off,  this  bit  being  the  only  difference  between  an  8-bit  ASCII 
character  and  a  UNIVAC  1100/82  9-bit  character. 

2.  Plot  Program  for  Photo  Data 

A.  Introduction.  This  section  describes  an  ASCII  FORTRAN  UNIVAC  1100/83® 
program  which  plots  UN IV AC-compatible  photo  data  as  (a)  X-Y  film  frame  plots 
and  (b)  X-Y  trajectory  plots.  These  plots  provide  a  method  for  error  detec¬ 
tion  of  faulty  photo  data,  allowing  for  correction  of  the  data  before  any 
further  processing  takes  place. 

B.  Discussion.  The  plot  program  generates  both  x-y  film  frame  plots  and  x-y 
trajectory  plots.  The  film  frame  plots  (see  Figure  5a)  plot  the  position  of 
each  photo  target  (on  a  reel  of  film)  relative  to  the  film  frame  origin.  Each 
individual  frame  of  film  is  plotted,  and  all  target  positions  in  the  frame  are 
connected  by  straight  lines;  this  aids  in  visual  recognition  of  relative 
changes  in  target  positions  from  frame  to  frame.  There  is  sometimes  more  than 
one  anatomical  mount  (see  Figure  6  for  mount  example)  visible  to  a  camera 
during  a  run.  It  is  for  the  reason  that  a  separate  set  of  film  frame  plots 
are  generated  for  each  anatomical  mount  found  on  a  reel  of  film.  For 
instance,  if  two  mounts  are  visible  to  a  particular  camera  and  only  one  mount 
Is  visible  to  a  second  camera,  three  sets  of  film  frame  plots  are  generated 
for  the  experiment.  A  total  of  250  film  frames  are  plotted  for  each  set. 


which  consists  of  ten  pages  of  twenty-five  plots  each.  On  top  of  each  page 
the  plot  scale  (meters)  is  listed;  this  measures  the  greatest  target  displace¬ 
ment  from  the  first  to  last  frame  on  the  page.  In  addition  all  target  posi¬ 
tions  which  were  manually  digitized  (not  automatically  tracked)  by  the  PDS 
operator  are  circled  for  identification.  The  film  frame  plots  are  used  to 
identify  errors  which  may  exist  in  the  data.  For  instance,  a  target  nay  have 
been  identified  incorrectly  by  the  PDS  operator;  this  problem  is  visually 
exposed  in  the  film  frame  plots.  Subsequently  the  data  file  is  corrected. 

The  x-y  trajectory  plots  (see  Figure  5b)  identify  the  path  each  target  (on 
a  reel  of  film)  travels  in  the  film  plane  during  a  run.  Each  individual  point 
on  the  plot  measures  target  displacement  from  a  fixed  film  frame  origin.  A 
total  of  250  data  points  are  plotted,  with  the  first  coinciding  with  time  of 
first  motion  of  the  experimental  subject.  Since  more  than  one  anatomical 
mount  may  appear  on  a  reel  of  film,  a  separate  x-y  trajectory  plot  is 
generated  for  each  mount  on  the  film.  The  path  of  every  photo  target  on  a 
mount  (visible  to  the  camera)  is  plotted,  each  identified  by  its  own  target 
number.  The  x-y  trajectory  plot  provides  at  a  glance  a  visual  presentation  of 
target  movement  (for  a  given  mount)  for  the  entire  run  duration;  this  allows 
for  quick  detection  of  data  discontinuities.  For  instance,  if  the  pattern 
recognition  algorithms  (used  in  the  photo  digitizing  process)  found  the  wrong 
target,  a  detectable  discontinuity  results  in  the  x-y  trajectory  plot. 

The  plot  program  for  photo  data  possesses  versatile  plot  output  capabilities, 
made  available  through  the  use  of  DISSPLA®  plot  software  (a  product  of 
ISSCO*).  All  plot  commands  are  output  to  a  catalogued  UNIVAC  disk  file  (via 
DISSPLA®  software).  This  "compressed"  plot  file  can  then  be  used  as  input  to 
a  variety  of  DISSPLA®  post-processor  programs,  which  "decompress"  the  plot 
commands  and  direct  them  to  a  plotting  device.  This  disk  file  may  be  utilized 
several  times  on  a  variety  of  plotting  devices  without  once  repeating  execu¬ 
tion  of  the  actual  plot  program.  ,  DISSPLA®  post-processor  programs  exist  for 
such  plotting  devices  as  the  TEKTRONIX  4014®,  FR80®  and  HP  7220®  (all  avail¬ 
able  to  NBDL) . 

There  are  advantages  in  using  any  of  these  post-processors.  For  instance,  the 
FR80®  post-processor  is  used  in  a  production  environment  to  generate  micro¬ 
fiche  and  high  quality  hardcopy.  The  TEKTRONIX  4014®  post-processor  can  be 
used  to  view  selected  plots  on  the  TEKTRONIX  4014®  scope  immediately  after 
execution  of  the  photo  data  plot  program.  This  quick  turnaround  is  of  benefit 
in  expedient  data  analysis.  The  FR80®  is  not  located  on-site,  so  in  order  to 
combine  the  benefits  of  quick  turnaround  and  high  quality  hardcopy  an  HP  7220® 
plotter  is  to  be  installed  on-site.  This  plotter  in  addition  affords  a  multi¬ 
color  plotting  capability. 


* 
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The  Input  to  the  plot  program  Is  UNIVAC-compatible  photo  data.  Character  data 
is  In  9-blt  UNIVAC  ASCII  format  and  numeric  data  Is  either  In  36-bit  UNIVAC 
Integer  format  or  36-bit  UNIVAC  floating  point  format  (depending  on  format 
specifications).  As  noted  In  section  l.B  data  Is  organized  on  the  Input 
files  by  film  reel,  with  a  separate  Input  file  for  each  camera  site.  A  direc¬ 
tory  exists  In  each  Input  file,  containing  the  run  numbers  (each  run  number 
Identifies  a  unique  film  reel)  and  accompanying  sector  locations  of  the  header 
and  time  records.  The  data  records  for  a  run  (or  film  reel)  Immediately 
follow  the  time  record,  and  a  "locator  table"  for  the  directory  precedes  all 
other  records  on  the  file.  Input  is  accomplished  with  a  general  purpose  I/O 
package  (see  Reference  1)  which  reads  FIELDATA  (NBDL  archival  standard)  In  an 
ASCII  FORTRAN  programming  environment.  Plots  are  output  for  each  Input  file 
(or  camera  site).  Within  each  file,  plots  are  output  by  run  number  (or  film 
reel)  and  within  each  run  plots  are  output  by  anatomical  mount. 

The  PDS  operator  is  responsible  for  entering  tlme-of-day  values  (LED  displays) 
for  each  frame  of  digitized  film.  To  insure  their  legitimacy  the  plot  program 
performs  tolerance  checks  for  the  time  values.  The  results  (see  Figure  5c), 
called  a  time  record  analysis,  are  output  to  a  line  printer*. 

The  plot  program  for  photo  data  accomplishes  the  following: 

(1)  Provides  error  detection  methods  for  photo  data  by  generating: 

(a)  x-y  film  frame  plots. 

(b)  x-y  trajectory  plots. 

(2)  Performs  confidence  checks  for  tlme-of-day  values  for  film  frames  and 
prints  the  results. 

(3)  Provides  plot  flexibility  features  through  the  use  of  post-processor 
programs. 

(a)  Plots  can  be  generated  repetitively  without  repeat  execution  of 
the  plot  program. 

(b)  Plots  can  be  viewed  Immediately  upon  execution  on  a  demand 
terminal  scope. 

(c)  Quality  hardcopy  and  microfiche  of  plots  are  generated. 

(4)  Accomplishes  Input  of  data  with  standardized  I/O  access  methods. 


*  This  analysis  may  be  output  to  a  plotting  device,  but  It  is  slower  and 
more  expensive  to  do  so. 


C.  Main  Program  and  Subroutine  Descriptions: 

(1)  PLOT  (see  Appendix  2b  for  listing)  -  This  Is  the  driver  program 
which  controls  all  program  operations.  A  single  Input  card  Identifies  (a)  the 
number  of  Input  files,  (b)  the  time  record  analysis  option,  (c)  the  X-Y  con¬ 
tour  plot  option,  (d)  the  frame  plot  option,  (e)  the  frame  plot  format  option, 
(f)  the  lowest  run  number  to  be  processed  and  (g)  the  highest  run  number  to  be 
processed  (see  Appendix  2a  for  card  format).  Once  the  Input  card  is  read, 
DISSPLA*  Is  called  to  initialize  the  compressor.  Plotting  Is  accomplished  by 
virtue  of  a  loop  which  controls  the  Input  and  subsequent  plotting  of  photo 
data  for  each  combination  of  run/camera  site/anatomical  mount  (see  Figures 
2b-2f  for  Input  record  formats). 

(2)  TIMER  (see  Appendix  2c  for  listing)  -  This  subroutine  generates  a 
time  record  analysis  for  a  given  run  number  and  site  ID.  The  time  record  Is 
Input  via  the  photo  data  I/O  program  (see  reference  1).  A  time  Interval 
between  frames  Is  approximated  by  averaging  total  elapsed  time  over  the  250 
frames  of  data.  This  Increment  Is  used  to  calculate  an  "expected  time"  for 
each  of  250  frames.  These  times  are  compared  to  the  actual  time  record  data 
points  In  order  to  test  for  excessive  error  (a  difference  between  actual  and 
"expected"  time  greater  than  .001  second).  Any  such  discrepancies  are  anno¬ 
tated  with  an  asterisk  (*)  on  the  time  record  analysis  (see  Figure  5c). 

(3)  LDATA  (see  Appendix  2d  for  listing)  -  This  subroutine  loads  photo 
displacement  data  Into  common  for  access  by  plot  subroutines.  This  routine  Is 
called  once  for  each  combination  of  run/camera/mount.  The  photo  data  I/O 
package  Is  used  to  accomplish  input  (in  order  to  read  FIELDATA  within  an  ASCII 
FORTRAN  program). 

(4)  MIXY  (see  Appendix  2e  for  listing)  -  This  subroutine  determines 
minimum  and  maximum  X  and  Y  valuqs  for  a  page  (twenty-five  frames)  of  frame 
plots.  These  values  are  placed  In  common  for  scaling  use  by  the  plot  subrou¬ 
tines.  The  minimum  and  maximum  values  are  selected  so  that  the  abscissa  and 
ordinate  axes  are  equal  lengths  In  data  units  (meters).  This  allows  the 
viewer  of  the  plots  to  determine  at  a  glance  In  which  direction  the  targets 
were  moving  fastest.  For  example.  If  the  targets  moved  from  the  far  left  of 
the  Initial  frame  on  the  page  to  the  far  right  of  the  final  (25th)  frame,  the 
X  fllmplane  displacement  was  changing  more  rapidly  than  the  Y  displacement. 

(5)  PLOTD  (see  Appendix  2f  for  listing)  -  This  subroutine  consists  of 
two  entry  points.  All  plotting  In  both  entry  points  Is  done  with  DISSPLA* 
software. 

A.  PL0T1  -  This  entry  point  plots  the  title  Information  for  each 
page  of  frame  plots.  Included  In  the  title  Information  am  the  (1)  run 
number,  (11)  camera  number,  (111)  plot  scale  (meters),  (1v)  anatomical  mount 
type  and  (v)  page  number.  All  plotting  Is  done  with  DISSPLA*  software. 


B.  PL0T2  -  This  entry  point  plots  all  frame  plots.  It  is  called 
once  for  every  frame  of  photo  data  to  be  plotted.  Two-hundred  fifty  frames 
are  plotted  for  each  mount,  twenty-five  frames/page,  ten  pages/mount.  The 
target  number  for  all  targets  tracked  (by  the  PDS  operator)  in  a  particular 
frame  are  printed  directly  at  their  digitized  position.  Those  targets  digi¬ 
tized  manually  are  circled  for  identification.  All  targets  in  a  particular 
frame  are  connected  by  vectors  in  order  to  present  the  image  of  motion  as  the 
targets  begin  to  change  position  from  frame  to  frame.  Targets  may  have  been 
deleted  or  added  as  the  frames  advance.  Deleted  target  numbers  are  printed  in 
the  bottom  right  corner  of  the  frame  in  which  the  target  was  dropped.  Added 
targets  are  printed  in  the  bottom  left  corner  of  the  frame  in  which  the  target 
was  added.  If  chair  data  are  available  for  the  current  run/camera,  a  * C ’  is 
printed  in  the  top  left  corner  of  each  frame.  The  frame  number  is  printed  in 
the  top  right  corner  of  the  frame. 

(6)  XYPLOT  (see  Appendix  2g  for  listing)  -  This  subroutine  plots  X  vs. 
Y  filmplane  displacement  values  for  a  given  run/camera/anatomical  mount.  All 
targets  digitized  for  the  mount  are  plotted  on  the  same  graph.  A  maximum  of 
250  data  points  are  plotted  (and  connected)  for  each  target,  with  the  target 
number  appearing  coincident  with  the  first  point  for  a  curve.  The  maximum  and 
minimum  values  of  the  abscissa  and  ordinate  are  printed  on  the  plot.  In  addi¬ 
tion,  titling  Information  is  listed  which  includes  (i)  run  number,  (1i)  camera 
number,  (iii)  anatomical  mount  type  and  (iv)  current  date.  The  X-Y  contour 
plot  affords  the  viewer  a  quick  method  of  isolating  data  incongruities. 

3.  Reformatting  (sign  change)  Program  for  Photo  Data 

A.  Introduction.  This  section  describes  an  ASCII  FORTRAN  UNIVAC  1100/83* 
program  which  converts  UN IV AC-compatible  photo  data,  digitized  manually,  into 
the  format  of  data  tracked  under  computer  control.  This  prepares  the  data  for 
subsequent  computation  of  3-dimensional  trajectories  of  test  subjects. 

B.  Discussion.  This  program  is  designed  to  read  all  target  position  data 
for  a  given  range  of  run  numbers  and  scan  for  negative  data  values.  Each 
negative  value  represents  an  entry  made  manually  by  the  PDS  (Photo  Digitizing 
System)  operator  using  the  x-y  crosshairs.  All  such  values  are  multiplied  by 
-1  and  rewritten  to  the  data  file,  ton  archival  data  tape  is  later  generated 
(with  the  reformatted  data  records)  by  an  independent  runstream.  This  tape 
serves  as  input  to  the  EASYFLOW  runstream. 

Input  to  the  reformatting  program  is  UNIYAC-compatible  photo  data.  Character 
data  are  in  9-bit  UNIVAC  ASCII  format  and  numeric  data  are  either  in  36-bit 
UNIVAC  Integer  format  or  36-bit  UNIVAC  floating  point  format  (depending  on 
format  specifications).  Data  are  organized  on  the  input  files  by  film  reel 
(or  run  number),  with  a  separate  Input  file  for  each  camera  site  (see  section 
l.B  for  a  more  detailed  description).  Direct  access  tables  accompany  each 
Input  file.  Input/output  is  accomplished  with  a  general  purpose  I/O  package 
(reference  1)  which  reads/writes  FIELDATA  (NBOL  archival  standard)  within  the 
context  of  an  ASCII  FORTRAN  programming  environment.  Output  file  formats  for 
the  reformatting  program  are  identical  to  the  Input  file  formats,  the  only 
difference  being  that  all  values  in  the  output  data  records  are  positive. 
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The  reformatting  program  accomplishes  the  following: 

(1)  Prepares  UNIVAC-compatlble  photo  data  for  subsequent  computation  of 
three-dimensional  trajectories  (of  test  sifcjects)  by  converting  manually 
digitized  photo  data  Into  the  format  of  computer-tracked  photo  data. 

(2)  Accomplishes  Input/output  of  data  with  standardized  I/O  access 
methods. 


C.  Main  Program  Description: 

SGFLIP  (see  Appendix  3b  for  listing)  -  This  Is  the  driver  program  which 
controls  all  program  operations.  The  only  subroutines  utilized  are  those 
belonging  to  the  photo  data  ASCII  FORTRAN  I/O  program.  A  single  Input  card 
Identifies  (a)  the  number  of  Input  files,  (b)  the  lowest  run  number  to  be 
processed  and  (c)  the  highest  run  number  to  be  processed  (see  Appendix  3a  for 
card  format).  The  number  of  Input  files  and  the  run  number  range  (from  the 
Input  card)  are  used  to  Index  a  loop  which  systematically  (a)  reads  data 
records,  (b)  scans  them  for  negative  data  points  and  (c)  changes  negative 
values  to  positive  values.  Once  a  data  record  Is  completely  processed  It  Is 
rewritten  to  the  Input  file  unless  no  negative  values  were  found  (see  Figures 
2b-2f  for  Input/output  record  formats). 

4.  Operation  of  Photo  Data  Reduction  System 

A.  Introduction.  This  section  describes  the  actual  UNIYAC  1100/83* 
production  run stream  which  controls  data  reduction  processing  of  NBDL  photo 
data  prior  to  Its  Input  to  the  EASYFLOW  runstream. 

B.  Discussion.  As  previously  noted,  there  are  three  main  processing  steps  In 
the  photo  data  reduction  process;  conversion,  plotting  and  reformatting. 

Each  step  performs  an  Integral  function  Independent  of  the  other  steps.  This 
modular  organization  of  steps  provides  several  benefits: 

(1)  Any  step  can  be  performed  Independently  from  the  others;  l.e.,  each 
step  has  stand-alone  capability. 

(2)  Modifications  to  steps  can  be  effected  with: 

(a)  Complete  non-interference  with  other  steps. 

(b)  Expediency  In  that  all  necessary  changes  take  place  In  a  single 

module. 

(3)  Input/output  may  be  examined  at  any  point  between  steps. 

(4)  Additional  steps  may  be  added  to  the  overall  process  without  altering 
any  existing  steps,  thus  providing  a  desirable  environment  for  enhancements  to 
the  system. 


9 


11 


Each  step,  or  module,  shall  heretofore  be  referred  to  as  an  element.  Element 
Is  a  UNIVAC  term  describing  a  collection  of  computer  commands.  When  it  Is 
desired  to  execute  an  element,  or  step,  the  computer  Is  Instructed  to  carry 
out  all  commands  In  that  element.  The  production  runstream  described  In  this 
section  executes  a  collection  of  three  such  elements: 

(1)  XQTCVT  Is  the  conversion  element. 

(2)  XQTPLT  Is  the  plot  element. 

(3)  XQTSGN  is  the  reformatting  (sign  change)  element. 

The  first  element  to  be  executed  Is  the  conversion  element  XQTCVT  (see 
Appendix  la).  The  main  function  of  this  element  Is  to  execute  the  photo 
conversion  program.  Additionally,  the  output  of  the  conversion  program  Is 
written  to  a  UNIVAC  data  tape.  In  sequence,  the  conversion  element  (1)  copies 
a  PDS  data  tape  onto  disk,  (2)  executes  the  photo  conversion  program  and  (3) 
writes  the  output  (from  the  conversion  program)  to  a  UNIVAC  data  tape.  The 
UNIVAC  tape  Is  saved  for  10  days  as  a  precautionary  measure.  The  key  output 
from  the  conversion  element  Is  the  converted  photo  data  written  to  disk  by  the 
conversion  program. 

If  It  Is  necessary  to  correct  for  PDS  operator  error  (see  section  l.B.)  the 
conversion  element  Is  altered  to  Include  correction  options.  These  options 
are  Input  to  the  element  In  card  Image  via  the  UNIVAC  EDIT  processor  (see 
reference  5).  The  EDIT  processor  commands  are  part  of  the  production 
runstream.  Thus  the  runstream  actually  edits  the  conversion  element  before 
executing  It.  The  correction  options  allow  for  rectification  (by  the 
conversion  program)  of  Instances  where  the  PDS  operator: 

(1)  Enters  Incorrect  run  numbers  on  PDS  header  records. 

(2)  Outputs  data  to  a  PDS  tape  for  two  film  reels  with  the  same  run 
number  and  camera  site  ID  (on  the  header  record). 

(3)  Enters  Incorrect  camera  site  ID's  on  PDS  header  records. 

The  first  card  Image  of  options  (Input  with  the  EDIT  processor)  Indicates  the 
number  of  film  reels  to  be  corrected  for  each  of  the  three  categories  of  pro¬ 
blem  cases.  A  card  Image  then  follows  for  each  occurrence  of  a  problem  case 
(1);  It  specifies  the  run  number  and  camera  site  ID  of  (the  header  record  for) 

the  problem  film  reel  as  well  as  the  occurrence  (of  this  unique  combination) 

to  be  corrected.  It  also  specifies  the  correct  run  number  to  replace  the 
Incorrect  field  In  the  header  record.  After  the  last  correction  card  image 
for  case  (1),  card  Images  follow  for  case  (2).  Each  of  these  specifies  the 

run  number  and  camera  site  ID  of  the  problem  reel.  The  first  occurrence  of 

this  combination  of  run  number  and  site  ID  (In  a  PDS  header  record)  Is  by¬ 
passed  and  only  the  second  occurrence  (of  PDS  records  with  this  header  Identi¬ 
fication)  Is  converted  to  UNIVAC  format.  Lastly,  the  correction  card  Images 
for  case  (3)  follow.  Each  specifies  the  run  number,  camera  site  ID  and 
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occurrence  Index  for  the  problem  reel.  Additionally,  It  contains  the  correct 
camera  site  10.  If  no  corrections  are  necessary  on  the  PDS  input  tape,  only 
one  card  image  is  input.  It  indicates  quantity  zero  for  each  of  the  three 
problem  cases. 

The  EDIT  processor  not  only  controls  input  of  all  correction  options  to  the 
conversion  element,  but  it  also  is  used  to  enter  the  reel  number  of  the  PDS 
data  tape  to  be  used  as  input  to  the  photo  conversion  program.  As  with  the 
EDIT  commands  for  correction  options,  these  commands  also  constitute  a  portion 
of  the  overall  production  runstream. 

The  second  element  of  the  runstream  to  be  executed  is  the  plot  element  XQTPLT. 
Its  primary  function  is  to  execute  the  plot  program  for  photo  data.  Addi¬ 
tionally,  the  disk  file  of  DISSPLA®  plot  commands  (output  from  the  plot  pro¬ 
gram)  is  used  as  input  to  the  FR80®  post-processor  program  to  produce  produc¬ 
tion  quality  microfiche  and  hardcopy  of  all  plots.  In  sequence,  the  plot 
element  (1)  executes  the  plot  program  for  photo  data  (the  plot  program  uses 
photo  data  residing  on  disk  as  Input),  (2)  executes  the  FR80®  post-processor 
program  and  (3)  writes  the  disk  file  of  plot  commands  to  a  catalogued  disk 
file  for  public  access.  As  previously  noted  the  plot  program  uses  output  from 
the  conversion  program  as  Its  input.  This  input  consists  of  disk  files  of 
photo  data,  one  for  each  camera  site  (1,  2  or  3)  used  in  a  series  of  runs. 
Output  from  the  plot  program,  a  disk  file  of  DISSPLA®  plot  commands,  serves  as 
input  to  the  FR80®  post-processor.  The  post-processor  interprets  the  plot 
commands  in  order  to  produce  microfiche  and  hardcopy  of  all  plots  generated  by 
the  plot  program  (these  include  the  x-y  film  frame  plots  and  x-y  trajectory 
plots).  These  plots  are  used  In  error  detection  analysis  and  are  maintained 
on-site. 

Certain  plot  options  exist  and  are  input  to  the  plot  element  via  the  UNIVAC 
EDIT  processor.  The  EDIT  commands  constitute  part  of  the  production  runstream 
(see  Figure  8);  they  are  used  to’  edit  the  plot  element  before  executing  it.  A 
single  card  Image  contains  all  the  plot  options.  These  options  Include: 

(1)  Number  of  Input  files. 

(2)  Time  record  analysis  option. 

(a)  Output  to  line  printer. 

(b)  Output  DISSPLA®  plot  commands  to  disk  file. 

(c)  Do  not  output  to  any  device. 

(3)  X-Y  trajectory  plot  option. 

(a)  Output  DISSPLA®  plot  commands  to  disk  file. 

(b)  Do  not  output  to  any  device. 
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(4)  X-Y  film  frame  plot  option. 


(a)  Output  DISSPLA*  commands  to  disk  file. 

(b)  Do  not  output  to  any  device. 

(5)  Plot  direction  option  for  x-y  film  frame  plots*. 

(a)  Plot  five  columns  to  a  page  with  frame  number  Increasing  In  each 
column  (plot  direction  always  downward). 

(b)  Plot  five  columns  to  a  page  with  frame  number  Increasing  In 
columns  1,  3  and  5  (plot  direction  downward)  and  decreasing  In  columns  2  and  4 
(plot  direction  upward). 

In  an  input  file  each  run 

(6)  Lowest  run  number  to  be  plotted  number  identifies  the  header 

record  for  a  single  film  reel; 

(7)  Highest  run  number  to  be  plotted  photo  data  Is  plotted  for  this 

range  of  run  numbers  (or  film 

reel s ) . 

Once  the  plot  program  has  been  executed,  the  plot  element  executes  the  FR80® 
post-processor  program  to  produce  microfiche  and  hardcopy  of  all  generated 
plots.  The  plot  (disk)  file  Is  then  copied  to  a  catalogued  disk  file. 

The  third  and  final  element  of  the  production  runstream  Is  the  reformatting 
element.  The  primary  function  of  this  element  Is  execution  of  the  refor¬ 
matting  (sign  change)  program.  The  element  also  generates  a  UNIVAC  data  tape 
containing  output  from  the  reformatting  program.  This  data  tape  Is  saved  for 
180  days  and  Is  used  as  Input  to  the  EASYFLOW  data  reduction  and  analysis 
system.  In  sequence,  the  reformatting  element  (1)  executes  the  reformatting 
program  and  (2)  saves  the  output  on  a  UNIVAC  data  tape.  The  reformatting  pro¬ 
gram  uses  photo  data  residing  on  disk  (output  from  the  conversion  program)  as 
input.  These  same  files,  after  being  reformatted  to  reflect  automatically 
digitized  photo  data,  are  output  to  the  UNIVAC  data  tape. 

The  UNIVAC  EDIT  processor  is  utilized  (In  the  production  runstream)  to  input 

(1)  reformatting  program  options  and  (2)  the  output  tape  file  name.  The  input 
options  include: 

(1)  Number  of  Input  files. 

In  an  Input  file  each  run 

(2)  Lowest  run  number  to  be  processed  number  identifies  the  header 

record  for  a  single  film  reel; 

(3)  Highest  run  number  to  be  processed  photo  data  is  reformatted  for 

this  range  of  run  numbers  (or 

film  reels). 


*  This  option  Is  used  only  In  the  event  x-y  film  frame  plots  are  generated. 


C.  Production  Runstream.  The  photo  data  reduction  runstream  Is  a  collection 
of  computer  commands  designed  to  execute  the  three  steps  (elements)  In  the 
data  reduction  process.  The  runstream  consists  of  (1)  UNIVAC  Operating  System 
commands  (see  reference  3)  and  (2)  UNIVAC  EDIT  processor  commands  (see  refer¬ 
ence  5).  All  three  elements  to  be  executed  are  stored  In  a  single  program 
file.  The  runstream  (1)  copies  the  program  file  onto  a  scratch  disk  file,  (2) 
enters  the  Input  options  for  each  element  and  (3)  executes  each  element.  The 
only  Input  to  the  runstream  Is  a  single  PDS  data  tape  containing  digitized 
photo  data.  Output  from  the  runstream  Includes: 

(1)  A  10-day  UNIVAC  data  tape  of  converted  photo  data. 

(2)  Catalogued  disk  file  containing  DISSPLA*  plot  commands  for  all 
generated  plots. 

(3)  Microfiche  and  hardcopy  of  all  generated  plots. 

(4)  A  180-day  UNIVAC  data  tape  of  converted  and  reformatted  photo  data. 

(5)  Line  printer  listing  containing: 

(a)  Time  record  analyses. 

(b)  Operating  system  diagnostics. 

(c)  User-generated  diagnostics. 

(d)  Operating  system  accounting  Information. 


Several  considerations  entered  Into  the  design  of  the  photo  data  production 
runstream.  The  more  significant  of  these  are  noted  as  follows: 

(1)  In  storing  plot  commands  on  a  catalogued  disk  file,  plots  can  be 
generated  as  many  times  as  necessary  (and  on  several  plotting  devices)  without 
once  repeating  execution  of  the  actual  plot  program. 

(2)  In  some  Instances  (a  film  reel  of)  photo  data  Is  redlgltlzed  In  order 
to  Improve  data  quality.  It  Is  Important  (In  computing  3-dlmenslonal  position 
of  the  test  subject)  that  a  minimum  total  of  three  targets  (not  all  from  the 
same  camera  site)  are  digitized  successfully  for  a  given  test  run;  these  tar¬ 
gets  may  be  tracked  from  any  of  the  available  film  reels  (camera  sites)  for 
the  run.  If  the  PDS  computer  cannot  successfully  locate  enough  targets 
(through  computer-controlled  tracking),  the  PDS  operator  generates  another 
data  tape  In  which  ( s) he  manually  digitizes  (with  the  x-y  crosshairs)  the 
necessary  targets. 


The  production  runstream  is  designed  to  reprocess  the  new  POS  data  tape  with  a 
minimum  of  operator  preparation.  It  Is  also  equipped  with  the  option  to 
process  redigitized  data  residing  on  the  same  PDS  tape  as  the  Initial  data 
(for  the  same  film  reel).  The  production  runstream  also  creates  UNIVAC  data 
tapes  directly  after  (a)  conversion  and  (b)  reformatting  in  order  to  make 
data  available  for  further  analysis,  error  correction,  and  reprocessing. 

(3)  The  individual  cameras  (used  In  the  experiments)  are  calibrated  to 
provide  mathematical  camera  constants  vital  In  the  computation  of  three- 
dimensional  position  (of  the  test  subject).  Examples  of  camera  constants  are 
the  camera  nodal  point  coordinates  (x,  y  and  z)  and  the  lens  distortion  para¬ 
meter.  These  values  and  other  camera  constants  are  generated  by  a  camera- 
calibration  program  which  uses  as  Input  converted  data  from  a  camera- 
calibration  film  reel.  The  photo  data  production  runstream  is  designed  to, 
with  slight  modification,  execute  the  camera-calibration  conversion  program. 
This  altered  runstream  is  Included  as  an  element  In  the  program  file  con¬ 
taining  all  photo  data  reduction  programs. 

(4)  The  modular  organization  of  the  steps  (elements)  In  the  photo  data 
reduction  process  allows  for  efficient  and  expedient  expansion  of  functions. 
New  elements  can  be  added  without  changing  any  existing  elements.  This 
feature  affords  the  system  a  desirable  flexibility.  Standardized  data  formats 
and  I/O  access  methods  provide  further  flexibility.  Possible  additional 
functional  elements  are  listed  below: 

(a)  Error  editing  elements  which  correct  improper  data  values  in  the 
data  file, 

(b)  Packing  elements  which  combine  data  files  (and  accompanying 
directories)  from  multiple  tapes  onto  a  single  tape  (with  one  all-inclusive 
directory) . 

(c)  Plot  elements  which  allow  plot  selection  of  individual  photo 
variables  residing  on  a  converted  photo  data  file. 

(d)  Special  analysis  elements  which  extract  the  data  easily  through 
the  use  of  standardized  data  file  formats  and  standardized  I/O  access 
methods. 

Operational  particulars  of  the  runstream  are  as  follows: 

( 1 )  Program  File: 

DATASYSTEMS*PHOTOSTREAM. 

(2)  Programming  Languages: 

(a)  UNIVAC  ASCII  FORTRAN  Level  I0R1 

(b)  UNIVAC  Executive  System  EXEC  Level  38R5 

(c)  UNIVAC  Text  Editor  Level  16R1 
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CONCLUSION 

The  photo  data  reduction  and  analysis  system  has  been  revised  to  successfully 
meet  all  proposed  objectives: 

1.  All  constituent  programs  are  now  compatible  with  announced  changes  In  the 
UNIVAC  operating  system.  These  changes  Incorporate  an  ASCII  FORTRAN  pro¬ 
gramming  environment. 

2.  A  new  graphics  support  capability  is  now  utilized  to: 

a.  Improve  error  detection  procedures  by  providing  additional  plot  pre¬ 
sentations. 

b.  Produce  output  compatible  with  a  variety  of  graphics  devices  (such  as 
the  Tektronix  4014®  and  the  Hewlett  Packard  7220®). 

c.  Reduce  turn-around  time  of  analysis  by  providing  for  Immediate  review 
of  plots. 

3.  A  stand-alone  runstream  providing  for  greater  operator  efficiency  is  now 
in  effect.  In  addition,  all  component  programs  exist  In  a  single  program 
file. 

4.  Standardized  data  formats  and  I/O  access  methods,  which  in  unison  provide 
compatibility  with  stored  archival  data,  are  now  incorporated  into  photo  data 
reduction  processing. 

5.  A  modular  organization  of  software  is  now  incorporated  in  order  to  provide 
for  greater  operational  flexibility. 


POS  OUTPUT  TAPE  FORMAT 


§ 

RECORD  SEQUENCE 


W 


Header 

Time 

Target  1, 

x  values 

N 

y  values 

Target  2, 

x  values 

II 

y  values 

Target  N, 

x  values 

M 

y  values 

} 

} 


End-of-FIle 


First  250 
frames 


Repeat  above 
format  if  more 
than  250  frames 
are  being 
digitized 


End-of-FIle 

End-of-FIle  (double  EOF:end  of  data  on  tape) 


Complete 
Date  for 
One  Film 


Repeat 

for 

each 

film 


i 


I 


March  29,  1979 
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Figure  la 


PHOTO  X  (OR  Y)  DATA  RECORD  FORMAT 


REPEAT  FOR  TOTAL  OF  UNIVAC  OUTPUT  BUFFER 

250  DATA  VALUES  U36)  ARRAY  LOCATION 


HEADER  FORMAT  (PDS  DATA) 


16  Bit  Word  Contents 

1  Record  type:  1  (16  bit  Integer) 

3>  Test  ID  (6  character  ASCII) 

4) 

5)  Time  of  first  frame  In  units  of  the  least  significant 

6 f  LED  digit  (32  bit  Integer) 

7  Julian  date  digitized:  DDDY  (16  bit  Integer) 

8  1000  X  machine  #  +  camera  site  #  (16  bit  integer) 

9  At:  expected  time  between  successive  frames.  In  units  of 

the  least  significant  LED  digit  (16  bit  Integer) 

10  Program  version  number  (16  bit  integer) 

11  Not  used 

12  Not  used 


March  29,  1979 


Figure  le 


TIME  RECORD  FORMAT  (PDS  DATA) 


16  Bit  Word 


Contents 

Record  type:  2  (16  bit  Integer) 

Test  ID  (6  character  ASCII) 

'M  TIME  MMM1  (12  character  ASCII) 

Time  of  frame  1  In  units  of  the 

least  significant  LED  digit  (32  bit  Integer) 


• 

509)  Time  of  frame  250 

510/ 

512  Not  used 

513  Not  used  (fixed  length  record) 


March  29,  1979 


Figure  If 
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X  (OR  Y)  COORDINATE  RECORD  FORMAT  (PDS  DATA) 


16  Bit  Word 


II 


5 

6 


9 

10. 


509) 

510/ 

511 

512 


Contents 

Record  type:  3*x  values  (16  bit  Integer) 
4*y  values 

Test  ID  (6  character  ASCII) 

Target  #  (32  bit  Integer) 


Julian  date  on  film:  DDDY  (16  bit  Integer) 
x  (or  y)  value  for  target  In  frame  1:  displacement  from  a 
fixed  location  In  the  fllmplane  in  units  of  .00001  Inches. 
Negative  value  Indicates  manual  entry  was  made  using  the 
crosshairs.  Zero  Indicates  location  of  target  was  not 
defined  for  this  frame.  (32  bit  integer) 


x  (or  y)  value  for  target  In  frame  250 
Not  used 

Not  used  (fixed  length  record) 


March  29,  1979 


Figure  lg 
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FILE  FORMAT  FOR  UNIVAC 


COMPATIBLE  PHOTO  DATA 


RECORD  SEQUENCE 


BOT 


Location  of  random  1 
Access  table  / 

Header  ) 

Time  J 

Target  1,  x  values  ' 
"  y  val ues 

Target  2,  x  values 
"  y  values 


5  words 
7  words 


253  words 
(each  record) 


Target  N,  x  values  I 
"  y  values  ) 

End-of-flle 

1  Repeat  above 
f  format  If  more 
)  than  250  words 
f  were  digitized 


} 

End-of-flle 

End-of-flle  > 


End-of-flle 

End-of-flle 

End-of-flle  (Triple  EOF:  end  of  data)  - 
Random  access  table 


Compl ete 

data 

for 

one 

run 


Repeat 

for 

each 

run 


SITE* 

ONE 


*  Above  block  of  data  copied  from  mass  storage  to  tape  using  @  copy. 

There  are  3  blocks  of  site  one  data  copied  to  tape  followed  by  3 
blocks  of  site  2  followed  by  3  blocks  of  site  3. 

May  5,  1982  Figure  2a 
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RECORD  FORMAT  FOR  LOCATION  OF  RANDOM 
ACCESS  TABLE  FOR  UNIVAC  COMPATIBLE  PHOTO  DATA 

WORD  CONTENTS 

1  Sector  address  of  random  access  table  (Integer) 

2  Number  of  runs  In  this  block  of  data  (Integer) 

3  Site  ID  (Integer) 

4  Not  used 

5  Not  used 


May  5,  1982 


Figure  2b 


HEADER  RECORD  FORMAT  FOR  UNIVAC 
COMPATIBLE  PHOTO  DATA 
CONTENTS 

Run  number  (alpha  -  6  characters)  FLOAT A 
The  word  'PHDATA*  (alpha)  FLDATA 
Site  ID  (integer) 

Time  of  first  frame  -  5  most  significant 
digits  (real) 

Julian  date:  DDDY  (alpha-6  char)  Date  dlgltf/ed 
At:  expected  time  between  frames  (real) 

Program  version  number  (alpha-6  char) 


Figure  2c 


TIME  RECORD  FORMAT  FOR  UNIVAC 
COMPATIBLE  PHOTO  DATA 
CONTENTS 

Run  number  (alpha  -  6  characters)  FLDATA 
The  word  'CTSbbb'  (alpha)  FLDATA* 

Sfte  ID  (Integer) 

Time  of  first  frame,  5  least  significant 
digits  (real) 


Time  of  250th  frame 


C  Is  the  camera  site  ID  (1,  2,  or  3) 


Figure  2d 


r! 

tir 

s 

►  * 

• 

r  ' 

r 

h ' 

DATA  RECORD  FORMAT  FOR  UNIVAC 

COMPATIBLE  PHOTO  DATA 

• 

WORD 

CONTENTS 

• 

P 

1 

Run  number  (alpha  -  6  characters) 

• 

2 

Displacement  ID  (alpha  -  6  characters) 

(see  Figure  4) 

1 

3 

Site  ID  (Integer) 

• 

i 

4 

• 

X  (or  Y)  value  for  target  In  frame  1:  Displacement  from  a 
fixed  location  In  meters.  Negative  value  Indicates  manual 
entry  was  made  using  the  crosshairs.  A  value  of  999.0 
indicates  location  of  target  was  undefined  for  this  frame. 

(REAL) 

•  ! 

a 

• 

• 

253 

X  (or  Y)  value  for  target  In  frame  250. 

p  ’ 

■ 

- —  -  •< 

9  . 

8. 

9  ... 

ft! 

•_  _ 

9l _ _ 

May  5,  1982  Figure  2e 
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RANDOM  ACCESS  DIRECTORY  FORMAT 


FOR  UNIVAC  COMPATIBLE  PHOTO  DATA 
FIELD  CONTENTS 

0  Words  following  In  this  record  (Integer) 

1  RN{1)  Run  number  for  first  run  (alpha  - 

6  characters)  FLDATA 

2  SH(1)  Sector  address  of  header  for  first  run 

(Integer) 

3  NR(1)  Number  of  physical  records  of  253  words 

required  to  make  one  logical  record  for  each 
target  for  first  run.  This  occurs  when  more 
than  250  frames  are  digitized  for  a  run.  Out¬ 
put  Is  written  In  multiples  of  253  word  physi¬ 
cal  records  (Integer) 

4  ST(1)  Sector  address  of  time  record  of  first 

physical  record  for  first  run  (Integer) 


NR(l)+3  ST(1,NR(1))  Sector  address  of  time  record  of  NR 

physical  records  for  first  run  (Integer) 


For  J  number  of  runs  where  X  Is  the  number  of  words  used  by  all 
previous  runs 

1+X  RN(J)  Run  number  for  Jth  run  FLDATA 

2+X  SH(J)  Sector  address  of  header  for  Jth  run 

3+X  NR(J)  Number  of  physical  records  for  Jth  run 

4+X  ST(0,1)  Sector  address  of  time  record  of  1st 

physical  record  for  Jth  run 


NR(J)+X+3  ST ( J ,NR( J ) )  Sector  address  of  time  record  of 

NR(J)th  physical  record  for  Jth  run 


See  Figure  3  for  example 


May  5,  1982 


Figure  2f 
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EXAMPLE  OF  RANDOM  ACCESS  TABLE  WITH 


SINGLE  PHYSICAL  RECORD  FOR  EACH  TARGET 


RUN 


NO. 

SH 

NR 

ST(1) 

LX3036 

1 

1 

2 

LX3037 

134 

1 

135 

LX3038 

267 

1 

268 

LX3039 

400 

1 

401 

LX3041 

533 

1 

534 

LX3042 

666 

1 

667 

LX3045 

799 

1 

800 

LX3047 

932 

1 

933 

LX3048 

1065 

1 

1066 

LX3050 

1198 

1 

1199 

LX3051 

1331 

1 

1332 

LX3053 

1504 

1 

1505 

LX3049 

1637 

1 

1638 

Figure  3 


p  I 


DISPLACEMENT  ID  (EZFLOW  FORMAT) 


CDMTTb 


6  char  ID  (Field  data) 


where 


camera  site  ID  (1,  2,  or  3) 
coordinate  component  (X  or  Y)  of  data 
mount  system 
1 :  neck* 


2:  mouth 


3:  top  of  head 
4:  pelvis 


5 :  other 


9:  sprocket  hole** 

TT:  target  ID  (01,  02,  03,  ...) 


*  The  chair  reference  target  has  the  special  mount/target  ID  of  100 

**  The  sprocket  hole  data  has  the  special  mount/ target  ID  of  900 


Figure  4 
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*TPF$(0). XQTCVT  (  26 ) 

•ASG.Td  1..U9H.REELN 


INPUT  DATA  FROM  DIGITI2ER 


•ASG  9. , F///300 

•ASG  10..F///3O0 

•ASG  II., F///300 

*. 

•XQT  MP0051/PHOTOC 

12345678901234567890123456789012345678901234567890 

». 


10 

•  .  1ST  CONTROL  STATEMENT 

1 1 

0. 

POS  5  - 

NRD  =  NO.  OF  RUNS  FOR  WHICH  TO  CHANGE  RUN  NOS.  MAX 

OF  6 

12 

P. 

POS  10  - 

NRB  =  NO.  OF  RUNS  TO  BYPASS  ON  THIS  TAPE.  MAX  OF  6 

13 

0. 

(SEE  NOTE  IN  PGM  MAIN) 

14 

P. 

POS  15  - 

NRC  =  NO.  OF  RUNS  REQUIRING  A  CHANGE  IN  CAM  NO.  MAX  OF 

15 

0. 

16 

p.  2ND  CONTROL  STATEMENT  (1  REQUIRED  FOR  EACH  MAGNITUDE  OF  NRD ) 

17 

0. 

POS  5-10 

-  LR  -  OLD  RUN  NUMBER 

18 

0. 

POS  15-20 

-  NR  *  NEW  RUN  NUMBER 

19 

0 . 

POS  25 

-  IC  =  CAMERA  NO.  FOR  WHICH  RUN  NO.  IS  TO  BE  CHANGEO 

20 

0. 

POS  30 

-  IP  *  OCCURRENCE  OF  RUN  NO.  TO  BE  CHANGED.  1ST. 

2ND. 

21 

0. 

MAY  BE  MORE 

THAN  ONE  OF  ABOVE  CONTROL  STATEMENT 

22 

0. 

23 

P.  NEXT  CONTROL 

STATEMENT  (1  REQUIRED  FOR  EACH  MAGNITUDE  OF  NRB) 

24 

0 . 

POS  5-10 

-  IRN  =  RUN  NUMBER  TO  BYPASS 

25 

0. 

POS  15 

-  ICN  =  CAMERA  NO.  FOR  WHICH  TO  BYPASS  RUN 

26 

0. 

MAY  BE  MORE 

THAN  ONE  OF  ABOVE  CONTROL  STATEMENT 

27 

0. 

28 

P.  NEXT  CONTROL 

STATEMENT  (1  REQUIRED  FOR  EACH  MAGNITUDE  OF  NRC) 

29 

p . 

POS  5-10 

-  LRN  =  RUN  NUMBER 

30 

p . 

POS  15 

-  LCN  «  OLD  CAMERA  NUMBER 

31 

P. 

POS  20 

-  NCN  =  NEW  CAMERA  NUMBER 

32 

P. 

POS  25 

-  NOC  =  OCCURRENCE  OF  CAMERA  TO  CHANGE,  1ST.  2ND, 

ETC. 

33 

P. 

MAY  BE  MORE 

THAN  ONE  OF  ABOVE  CONTROL  STATEMENT 

34 

P. 

35 

•  FREE 

1 . 

36 

•TEST 

TNE/0/S6 

37 

•JUMP 

OUT 

38 

•ASG.TF  A..U9S 

, SAVEOI  PREPLOTPHOTODATA 

39 

•COPY 

GM  9 .  ,  A  . 

40 

•TEST2 : 

41 

•TEST 

TG/1/S6 

42 

•JUMP 

TEST3 

43 

•COPY 

GM  10.,  A. 

44 

•TESTS: 

45 

•TEST 

TG/2/S6 

46 

•JUMP 

OUT 

47 

•COPY 

GM  1 1 .  , A . 

48 

•OUT: 

49 

•FREE 

A  . 

50 

•.  EXITED  FROM 

XQTCVT 

i  . 


•PRT.S  .MAIN 
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LAMBERT *TPF$(0). MAIN( 37) 

1  C  THIS  PROGRAM  WILL  READ  DIGITIZED  PHOTODATA  IN  PDS  OUTPUT 

2  C  FORMAT  AND  CONVERT  TO  UNIVAC  1182  F IELDATA  FORMAT.  A 

3  C  SEPARATE  FILE  IS  GENERATED  FOR  EACH  CAMERA  SITE  ID. 

4  C 

5  C  NRD  -  NUMBER  OF  RUNS  REQUIRING  A  CHANGE  IN  THE  RUN  NUMBER. 

6  C  THE  RUN  NUMBER  MAT  OR  MAY  NOT  BE  DUPLICATED. 

7  C  LR  =  OLD  RUN  NUMBER 

8  C  NR  *  NEW  RUN  NUMBER 

9  C  IC  *  CAMERA  NUMBER  FOR  WHICH  A  RUN  NUMBER  IS  TO  BE  CORRECTED 

10  C  IP  =  OCCURRENCE  OF  RUN  NUMBER  TO  BE  CHANGED.  1ST.  2ND.  ETC. 

1 1  C 

12  C 

13  C  NRB  =  NUMBER  OF  RUNS  TO  BYPASS  ON  THIS  TAPE 

14  C  IRN  =  RUN  NUMBER  TO  SKIP 

15  C  ICN  =  CAMERA  NUMBER  FOR  WHICH  TO  SKIP  RUN 

16  C 

17  C  IT  IS  ASSUMED  THAT  WHEN  A  RUN  IS  TO  BE  BYPASSED.  IT  HAS 

18  C  BEEN  PREVIOUSLY  DIGITIZED  AND  PLACED  IN  THE  FILE.  AND 

19  C  LATER  RE-DIGITIZED  AND  PLACED  IN  THE  FILE  AFTER  A  GROSS 

20  C  ERROR  WAS  DISCOVERED.  THE  FIRST  OCCURRENCE  OF  THE  RUN  IS 

21  C  ALWAYS  BYPASSED. 

22  C 

23  C 

24  C  NRC  «  NUMBER  OF  RUNS  REQUIRING  A  CHANGE  IN  CAMERA  NUMBER 

25  C  LRN  *  RUN  NUMBER 

26  C  LCN  =■  OLD  CAMERA  NUMBER 

27  C  NON  «  NEW  CAMERA  NUMBER 

28  C  NOC  =  OCCURRENCE  OF  CAMERA  NUMBER  TO  BE  CHANGED.  1ST.  2ND.  ETC 

29  C 

30  C 

31  DIMENSION  LCN(6),  NCN(6),  N0C(6),  ITRY(6),  ICCD(6) 

32  DIMENSION  IC(6).  IP(6).  ICN(6),  IPAS(6) 

33  DIMENSION  I TABLE ( 5 ) , dT ABLE ( 4000) . ICAM( 3 ) 

34  COMMON/ PHTHOR/RUN 10 , VAR 10 .ISITID.F1T1ME. dDTDIG, DELTA , IPGVER 

35  COMMON/ 1 36/ 1 OUT ( 3 ) . OATA( 250) . IDL . IUNIT , ITARG 

36  CHARACTER*6  RUN ID, VAR  ID , LR(6),NR(6),IRN(6), FD2ASC , 

37  1LRNC6) .dDTDIG. IPGVER. RUN< 100.3) 

38  C 

39  C 

40  DATA  IOUT/9. lO, 1 1/ 

41  68  FORMAT! IX. 14) 

42  I ERR*0 

43  I UN I T  *  1 

44  C 

45  C  READ  CORRECTION  INFORMATION 

46  C 

47  READ  (5.230)  NRD, NRB, NRC 

48  PRINT  230.  NRD. NRB. NRC 

49  C 

50  IF  (NRD.EQ.O)  GO  TO  20 

51  DO  10  d* 1 . NRD 

52  READ  (5,240)  LR( d ) , NR( d ) . IC ( d ) , I P( d ) 

53  PRINT  240,  LR( d ) . NR( d ) . IC(J ) . I P( d ) 

54  10  CONTINUE 

55  C 

56  20  IF  (NRB.EQ.O)  GO  TO  40 

57  DO  30  d*  1  .NRB 

58  RE AO  (5.250)  IRN( d ) . ICN( d) 

59  PRINT  250,  IRN( d ) . ICN( d ) 
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60 

30 

CONTINUE 

61 

C 

62 

40 

IF  (NRC.EO.O)  00  TO  60 

63 

00  50  J* 1 , NRC 

64 

READ  (5,250)  LRN( J) . LCN( d ) . NCN( J ) . NOC( J ) 

65 

PRINT  250,  LRN( J) , LCN( J ) ,NCN( J ) , NOC( J ) 

66 

50 

CONTINUE 

67 

C 

68 

60 

CONTINUE 

69 

C 

70 

NC*0 

71 

C 

72 

70 

CONTINUE 

73 

C 

74 

C 

READ  RECORO  FROM  POS  TAPE 

75 

C 

76 

CALL  CONVT ( I T YPE , I ST ) 

77 

IF  (IST.NE.4)  00  TO  210 

78 

C 

79 

IF  (ITYPE.NE.1)  GO  TO  70 

80 

IF  (NRD.EO.O)  GO  TO  100 

81 

DO  90  d-I.NRD 

82 

IF  (RUNIO . EQ. LR( d) . AND . ISITID. EQ.IC(d))  I TRY (d)«ITRY( J )+ 1 

83 

IF  ( RUNID . EO. LR( J ) . AND .ISITID.EQ.IC(J). AND . ITRY ( J) . EO ■ IP( J ) ) 

84 

1  RUNID'NR(J) 

85 

90 

CONTINUE 

86 

C 

87 

IOO 

IF  (NRB.EO.O)  GO  TO  120 

88 

DO  HO  d-I.NRB 

89 

IF  (RUNID. EQ.lRN(d). AND. ISITID.EQ. ICN( J) . AND . IPAS( J) . EO.O) 

90 

1  GO  TO  190 

91 

110 

CONTINUE 

92 

C 

93 

C 

PROCESS  RUNS  CONTAINING  INCORRECT  CAMERA  NUMBERS 

94 

C 

95 

120 

IF  (NRC.EO.O)  GO  TO  140 

96 

DO  130  J* 1 ,NRC 

97 

IF( RUNID. EQ.LRN(d).ANO. ISITID. EQ.LCN(J)) 

98 

1  ICCD(J)>ICCD(J)+1 

99 

IF  ( RUNIO. EQ.LRN(J). AND. ISITID . EQ. LCN(d ). AND. NOC( d).EQ. 

ioo 

1  ICCD(d) )  ISITID-NCN(d) 

101 

130 

CONTINUE 

102 

C 

103 

140 

CONTINUE 

104 

C 

105 

ISW'ISITIO 

106 

IF  (1SW.LE.0.0R.ISW.GT.3)  GO  TO  150 

107 

C 

108 

C 

DETERMINE  WHICH  CAMERA  NUMBER  TO  PROCESS 

109 

C 

110 

GO  TO  (160.170.180).  ISW 

111 

C 

117 

150 

CONTINUE 

113 

PRINT  270.  RUNIO. ISITID. ISW 

114 

STOP  'CAM  NO.  ERROR' 

115 

C 

116 

160 

CONTINUE 

117 

C 

118 

C 

CONVERT  DATA  FOR  CAMERA  1 

119 

C 

t 


9 


» 


r 
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CALL  CNVT 1(1  TYPE .1ST) 
IF  (NC.LT.2)  NC* 1 
00  TO  70 


i 


120 

121 

122 

123 

124 
129 
129 
127 
129 

129 

130 

131 

132 

133 

134 
139 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 
1  17 


c 

170  CONTINUE 
C 

C  CONVERT  DATA  FOR  CAMERA  2 
C 

CALL  CNVT  MITVPE,  1ST) 

IF  (NC.LT.3)  NC-2 
GO  TO  70 
C 

180  CONTINUE 
C 

C  CONVERT  DATA  FOR  CAMERA  3 

C 

CALL  CNVTKITYPE,  1ST) 

NC-3 

GO  TO  70 
C 

190  CONTINUE 

IPAS(d)»1 

C 

C  BYPASS  CURRENT  RUN 

C 

200  CALL  CONVT ( I TYPE , 1ST ) 

IF  (IST.E0.4)  GO  TO  200 
GO  TO  70 


148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 


C 

210  CONTINUE 
C 

C  WRITE  TRIPLE  END-OF  -FILE  AT  END  OF  DATA 

C 

DO  220. 1-1. NC 

OO  225. J-1.3 

CALL  EOFPUT ( I OUT ( I ) . IER ) 

IF(IER.NE.O) THEN 
PRINT  13. IER 
STOP  'ABORT  IN  MAIN’ 

END  IF 

229  CONTINUE 
220  CONTINUE 
C 

C  SET  CONDITION  WORO  TO  NUMBER  OF  CAMERAS 
C 

CALL  FCONO(ICONO) 

BITS! ICOND . 25 . 12) -BI TS( NC ,  25 , 12) 

CALL  FSETC(ICONO) 

C 

C  PRINT  DIRECT  ACCESS  TABLES 
C 

PRINT  251 
MAX-4000 

DO  300.I-1.NC 

CALL  GETDIR( IOUT( I ) .MAX ,0, ITABLE .IER) 

I F ( I ER . NE . O ) STOP  'DIRECTORY  ERR' 

ICAM( I )-ITABLE(3) 

NRUNS-ITA8LE(2 ) 

ISECT-ITABLE(I) 

PRINT  211, ICAM( I ) 


ISO 

CALL  GETOIR( IOUT( I ) .MAX, I SECT, JTABLE, IER) 

181 

IFOER.NE  .0)ST0P  'DIRECTORY  ERR' 

182 

IPRINT*2 

183 

DO  296 ,d*1, NRUNS 

184 

RUN( NRUNS . ICAM( I ) ) *F02ASC( JTABLE ( IPRINT ) , 6 ) 

185 

PRINT  22 1 , RUN( NRUNS . ICAM( I ) ) , ( JTABLE ( K ) , K- 1 PRINT  + 1 , IPRINT+1+ 

186 

1 dT ABLE (IPRINT+2)+1) 

187 

IPRINT“IPRINT+3+dTABLE( IPRINT+2) 

188 

296 

CONTINUE 

189 

PRINT  23 1 , ( ITABLE(K) ,  K»  1 ,5) 

190 

300 

CONTINUE 

191 

C 

192 

STOP 

193 

C 

194 

13 

FORMAT( IX, 'STATUS  ERROR  =  '.03) 

195 

230 

FORMAT  (315) 

196 

240 

FORMAT  ( 2( 4X , A6 ) , 215 ) 

197 

250 

FORMAT  (4X.A6.315) 

198 

260 

FORMAT  ('  READ  ERROR.  STATUS  =  '.13) 

199 

270 

FORMAT  ('  CAMERA  NUMBER  ERROR,  RUN  *  '.A6.2I10) 

200 

251 

FORMAT ( ' 1 ' ) 

201 

211 

FORMAT ( '  DIRECT  ACCESS  TABLE  FOR  CAMERA  '.12) 

202 

231 

FORMAT ( IX. S< 16. 4X)//) 

203 

221 

FORMAT ( 1X.A6.2X, 1318) 

204 

C 

205 

END 

•PRT . S 

. CONVT 
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37 
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TPF$(0 )  • CNVT 1(7) 

SUBROUTINE  CNVT 1 ( ITVPE . 1ST) 

C 

C  FUNCTION: 

C  CONTROL  (1)  READ  OF  DIGITIZED  PHOTO  DATA  FROM  PDS  TAPE 

C  (2)  CONVERSION  TO  UNIVAC  COMPATIBLE  FORMAT  AND  (3)  SCALING 

C  OF  DATA 

C 

C  ARGUMENT  DEFINITIONS: 

C  ITVPE -RETURNED  POS  RECORD  TYPE  WHERE 

C  1 “HEADER  RECORD 

C  2 “TIME  RECORD 

C  3=X  DATA  RECORD 

C  4=Y  DATA  RECORD 

C  IST-RETURNED  ERROR  STATUS  WORD  FROM  ADTIO  GENERAL 

C  PURPOSE  I/O  ROUTINE 

C 

COMMON/PHTHOR/RUNIO, VAR  ID, IS  IT  ID, FITIME, JDTOIG, DELTA, IPGVER 
COMMON/136 /I OUT (3) .DAT (250) . IDL , IUNIT , ITARG 
CHARACTER*6  VAR1C1(28)  .  VAR2C1(28),  VAR3C 1(28),  VARCC1(28) 
CHARACTER'S  VAR1C2(28),  VAR2C2(28),  VAR3C2(28).  VARCC2(28) 
CHARACTERS  VAR3C3(28) 

CHARACTER*6  RUNIO . VAR  10 . JDTOIG . IPGVER , B( 28 ) . TARNUM'3 , CAMNO* 1 
1 .OLDRUN 
C 

C  DEFINE  MOUNT/TARGET  ID'S 

C 

DATA  VAR  1C  1  / ' IX 100' . ' 1 Y 100' . ' IX 102 * . ' 1 Y 102 ' , 

1 ' 1X103' . ' 1 Y 103 ' .' 1X101 ',' 1Y 101 ', ' 1X201 ' , ' 1Y201 ' , 

2' 1X204' , ' 1Y204' , ' 1X203' , ' 1Y203' , ' 1X308' . ' 1Y30B' . 12*'  '/ 

DATA  VAR2C1  / ' IX lOO' , ' 1 Y lOO' . ' IX 103 ' , ' 1 Y 103 ' . ' IX 105 ' , 

1 ' 1Y105' , ' 1X104' , ' 1Y104' . ' 1X101' , ' 1Y101 ' , ' 1X201’ , ' 1Y201 ' , ' 1X204' 
2 ' 1 Y204 ' , ' 1X205' . ' 1Y205' . ' 1X203' . ' 1Y203' . ' 1X3 10' , ' 1 Y3 10' . 8* '  '/ 
DATA  VAR3C1  / ' IX 100' , ' 1 Y 100' . ' 1X2 12 ' . ' 1 Y2 12' . ' 1X204 ' , 

1 ' 1 Y 204 ' , '  1X208' . ' 1Y208' . ' 1X210' . ' 1Y210' , '1X202' , ' 1V202', '1X211' 
2' 1Y211 ' , ' 1X209' . ' 1 Y 209 ' . ' 1X207' , ' 1Y207' , ' 1X010' . ' 1 YOlO' , 8* '  '/ 
DATA  VARCC1  /' IX 100’ ,' 1Y 100' , ' 1X002 ',' 1Y002 ',' 1X003 ' , 

1 ' 1 Y003 ' , ' 1X004' , ' 1Y004' . ' 1X005' . ' 1Y005' . ' 1X201 ' , ' 1Y201 ' , ' 1X204' 
2 ' 1 Y204 ' , ' 1X205' , ' 1Y205' , ' 1X203' , ’ 1Y203' , ’ 1X310' , ' 1Y310' .8*'  '/ 
DATA  VAR1C2  /' 2X 100' ,' 2Y lOO' 2X203 ',' 2Y203 ' , 

1 ' 2X204 ' , ' 2Y204 ' , '2X202' . '2Y202' ,20* '  '/ 

DATA  VAR2C2  /' 2X 100 ',' 2Y lOO' ,' 2X205 ',' 2Y205 2X204 ’ , 

1 ' 2  Y 204 ' , '2X202' , '2Y202' . '2X203' . '2Y203' , 18* '  '/ 

DATA  VAR3C2  / ' 2X 1O0 ' . ' 2Y 100' . ' 2X 103 ' , ' 2Y 103 ' , 

1 '2X105' , '2Y105' , '2X101 ', '2Y 101' . '2X102' , 

2 ' 2Y 102 ' , '2X201 ' , '2Y201 ' . '2X202' , '2Y202' , '2X205' . '2Y205' , '2X204' 
3 ' 2Y204 ' , '2X203' , '2Y203' , '2X011' , '2Y011' , '2X012' , '2Y012' .4*'  '/ 
DATA  VARCC2  / ' 2X lOO ' , ' 2 Y lOO' . ' 2X 103 ' , '2Y 103 ' , ' 2X 105 ' , ' 2Y 105 ' , 

1 '2X101 ', '2Y 101 '. '2X102' , '2Y 102' . '2X201 ' . '2Y201 ' , '2X202' . 

2 ' 2Y202 ' , '2X205' . '2Y205' , '2X204' . '2Y204' , '2X203' . '2Y203' ,8* '  '/ 
DATA  VAR3C3  / ' 3X 100' . ' 3Y 100' , ' 3X 106 ' . ' 3Y 106 ' , ' 3X 1 1 1 ' , ' 3Y 1 1 1 ' , 

1 '3X102' , '3Y 102' . '3X101 ' . '3Y1Q1' , '3X103' . '3Y 103' . '3X107' , '3Y107' 
2 ' 3X 108 ' , ' 3Y 108 ' , '3X105' , '3Y 105' , '3X104' . ' 3Y 104 ' . '3X109' , '3Y109' 
3'3X1 10' . '3Y 1 10' ,4*'  '/ 

C 

OLDRUN* ' 

C 

C 

C  CONVERTS  AND  SCALES  PHOTO  DATA 
C 
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60 

6t 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


DO  90  KP>1 .  1000 
DO  lO  d- 1.250 
10  DAT ( 0 ) *0 . 0 

IF  (ITVPE.NE. 1)  GO  TO  20 
C 
C 

C  CHECK  FOR  MORE  THAN  250  FRAMES  OF  DATA. 

C 

C 

IF  (RUNID. EO.OLDRUN)  GO  TO  30 

OLORUN*RUNID 

VARID* 'PHOATA ' 

C 

C  WRITE  HEADER 

C 

CALL  PUTRCO( IOUT( ISITID) , RUNID , VARID. ISITID, IDL.DAT. IER) 

IF( IER.NE . 0)G0  TO  112 

PRINT  190, RUNID. VARID. ISITID. FITIME, JDTDIG, DELTA. IPGVER.ITYPE 
C 

C  READ  NEXT  PDS  RECORD 

C 

20  CALL  CONVT ( I T VPE . I ST ) 

IF  (IST.NE.4)  GO  TO  100 
I F ( I TYPE . EQ. 1 )G0  TO  90 
GO  TO  40 
C 
C 

C  HANDLES  EOF'S  FOR  MORE  THAN  250  FRAMES 

C 

30  CALL  EOFPUT( IOUT( ISITID), IER) 

IF ( IER. NE . 0)GO  TO  112 
CALL  CONVT ( I TYPE , 1ST ) 

IF  (IST.NE.4)  GO  TO  100 
C 
C 

C  DETERMINE  VARIABLE  NAMES  FOR  A  GIVEN  CAMERA. 

C 

C 

40  CONTINUE 

C  CHANGE  DATA  TO  METERS  8  SET  ZERO  VALUES  IN  NON-TIME  RECORDS  EQUAL 
VARID* ' 

ENCODE  ( 1 . 160.CAMN0)  ISITID 

SUBSTR( VARID. 1 , 1 )*CAMNO 

IF( I TYPE . EQ. 2)SUBSTR( VARIO, 2,2)*'T5' 

IF  (ITYPE.EQ.2)  GO  TO  70 
IF  (BITS( ITARG.5, 16) .NE .0)  GO  TO  60 
IF  (RUNID.LT. 'LX0253' )  NOT-1 
IF  (RUNID.LT. 'LX1410' )  NDT-2 
IF  (RUNID. GE. 'LX1410' )  NOT-3 

IF  (( RUNID. EQ. 'LX1234' ) .OR. (RUNID. EQ. ' LX  1235' ) .OR. (RUNID. EQ. 
1'LX1236' ). OR. (RUNID. EO. 'LX1237'))  NDT-4 
DO  50  d* 1 . 28 

IF  ((ISITID. EQ. 1). AND. (NOT. EQ. O)  B(d)*VAR1C1(d) 

IF  ((ISITID. EQ.1). AND. (NOT. EQ. 2))  B(d)-VAR2C1 (d) 

IF  ((ISITID.  EQ.D.ANO.  (NOT.  EQ.  3))  B(d)-VAR3C  1  ( d ) 

IF  ((ISITID. EQ.  1). AND. (NOT. EQ. 4))  B(d)-VARCCKd) 

IF  ((ISITI0.EQ.2) .AND. (NOT .EQ. 1 ) )  B(d)-VAR1C2(d) 

IF  (( ISITID. EQ. 2) .AND. (NOT. EQ.2) )  B(d)-VAR2C2(d) 

IF  ((ISITID. EQ.2). AND. (NOT. EQ. 3))  B(d)-VAR3C2(d) 

IF  ((ISITID. EQ.2). AND. (NOT. EQ. 4))  B(d)*VARCC2(d) 
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120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 


60 


70 


80 


IF  ((ISITID.EQ.3). AND. (NOT . EO. 3) )  B(d)«VAR3C3( d) 
CONTINUE 
ITT*2*ITARG 

IF  (ITYPE.E0.3)  VARID*B( ITT  -  1 ) 

IF  (ITYPE. EQ. 4)  VARID*B( ITT) 

GO  TO  70 
CONTINUE 

I F ( I TYPE . EO. 3)SUBSTR( VARIO , 2 , 1)«'X' 

IF( ITYPE . EQ. 4 )SUBSTR( VARID , 2 , 1)-'Y' 

NTC*BITS( ITARG, 5 , 16) 

ENCODE  (3, 1 70 , T ARNUM )  NTC 
SUBSTR( VARID , 3 . 3 ) =T ARNUM 
CONTINUE 
DO  80  d* 1.250 

IF  (ITYPE. EO. 2)  DAT(d ) *DAT ( d  )/ 10000 . 0 
IF  (ITYPE. GT. 2)  DAT(d ) =(DAT( d )/ 100000. 0)* . 0254 
IF  ((ITYPE. GT. 2). AND . (OAT (d).EQ.O.O))  DAT(d)-999.0 
CONTINUE 

IF  ( ( RUNID . GE . ' LX0253 ' ) . AND . ( VARID . EO . '  '))  GO  TO  90 


C 

C  CALL  SPROCKET  HOLE  SUBROUTINE  IF  900  OR  100  RECORD 
C 

IF  ( SUBSTR( VARID , 2,5). EO. 'X900  ')  CALL  SPH9X  (DAT) 

IF  (SUBSTR( VARID, 2.5) .EO. 'Y900  ')  CALL  SPH9Y  (DAT) 

IF  ( SUBSTR( VARID, 2 , 5 ) . EO- 'XtOO  ')  CALL  SPH1X  (DAT) 

IF  ( SUBSTR( VARIO , 2 . 5 ) . EO . ' Y 100  ')  CALL  SPH1Y  (DAT) 

C 

C 

C  WRITE  TIME  OR  DISPLACEMENT  RECORD 

C 

CALL  PUTRCD( IOUT(ISITIO). RUNID. VARID. ISITID, IDL. DAT. IER) 
IF ( IER . NE . 0)GO  TO  112 

PRINT  200. RUNID. VARID. ISITID. ITYPE. (DAT(d1 ).d1*1. 10) 

90  CONTINUE 
lOO  CONTINUE 


C 

C  WRITE  DOUBLE  ENO-OF-FILE  AT  END  OF  RUN  (END  OF 
C  DATA  FOR  1  FILM) 

C 

CALL  EOFPUT( IOUT( ISITID), IER) 

IF( IER . NE . O )G0  TO  112 

CALL  EOFPUT( IOUT( ISITID), IER) 

IF( IER . NE . 0)G0  TO  112 
C 

C  CALCULATE  SPROCKET  HOLE  AVERAGES 

C 

CALL  SPHAVG  (SPHX.SPHY) 

PRINT  180.  SPHX.SPHY 
C 

RETURN 


C 

C  OUTPUT  ABORT  ROUTINE 
C 

112  PRINT  12, IER 

12  FORMAT ( IX. 'STATUS  ERROR  «  '.03) 

STOP  'ABORT' 

C 

160  FORMAT  (II) 

170  FORMAT  (13) 

180  FORMAT  (/'  SPROCKET  HOLE  AVERAGES,  X  -  '.F10.6,'  Y  ■  '.F10.6) 


P 


47 


LAMBERT* TPF$(0 ) . CONVT (  34 ) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

1 1 
12 

13 

14 

15 

16 

17 

18 
19 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  CONVT ( I  TYPE , I  ST  ) 

FUNCTION: 

(1)  READ  POS  PHOTO  DATA  FROM  TAPE  AND  (2)  CONVERT 
OIGITIZEO  PHOTO  DATA  TO  UHNIVAC  COMPATIBLE  FORMAT 

ARGUMENT  DEFINITIONS: 

I TYPE -RETURNED  PDS  RECORD  TYPE  WHERE 
1 “HEADER  RECORD 
2*TIME  RECORD 
3=X  DATA  RECORD 
4  =  Y  DATA  RECORD 

IST-RETURNED  ERROR  STATUS  WORD  FROM  ADTIO  GENERAL 
PURPOSE  I/O  ROUTINE 

DIMENSION  132(300) 

COMMON/PHTHDR/RUNIO ,VAR ID, ISITID, FITIME, UDTOIG , DELTA , 1 PGVER 
COMMON/136/ IOUT( 3), DAT ( 250 ) , IDL , IUNIT , ITARG 
CHARACTERS  RUNID . VARIO , JDTDIG , I PGVER 


20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


C 

C  DATA  FOR  EXTENDING  THE  SIGN  BIT. 

C 

DATA  NZ  /07 40000000000/  .SIGN  /0020000000000/ 

C 

NW=255 

CALL  ATRD( IUNIT. 132. NW, 1ST) 

IF  (IST.EO. 1)  RETURN 
C 

C  ZERO  OUT  RECEIVING  ARRAY. 

C 

DO  10  1=1.250 
DAT ( I ) =0 . 
lO  CONTINUE 

C 

C  CONVERT  RUN  NUMBER  FROM  NOVA  ASCII  TO  UNIVAC  ASCII 
C 

CALL  ASC8T9( 132,1,17,6, RUNID ) 

C 

ITYPE=BITS( 1 32 (  1)  . 1 , 16) 

IF  (ITYPE.GT. 1)  GO  TO  60 
C 

C  EXTRACT  HEADER  RECORD  AND  CONVERT  TO  UNIVAC  WORDS 

C 

BITS(ITIME.5,8)=BITS(I32(2),29,8) 

BITSCITIME, 13,24 )*BITS( 132(3), 1.24) 

F 1  TIME “FLOAT ( ITIME ) 

BITS(ISITID,21, 16)=BITS(I32(4),5, 16) 

BITS( IDATE ,21,  1 2  )  =BI TS( 1 32 ( 3 ) , 25 , 1 2 ) 

BITS( I DATE ,33.4)*BITS(I32(4), 1,4) 

ENCODE (6.  151, JDTDIG) IDATE 

BITS(I DELTA, 33, 4)=BITS(I32(5), 1,4) 

OELTA=FLOAT(IOELTA) 

BITS(IVER,25. 1 2 ) =B I TS( 1 32 ( 5 ) ,5, 12) 

ENC0DE(6, 151 , IPGVER)IVER 
151  FORMAT ( 16) 

IF( IDATE. EQ.O) JDTDIG*' 

IF(IVER.EO.O)I PGVER* ' 

C 

RETURN 
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60 


IF  (ITYPE.GT.2)  GO  TO  70 
GO  TO  80 


LAMBERT*TPF$(0) .SPHAVG( 1 1 ) 

1  SUBROUTINE  SPHAV6 

2  DIMENSION  X(1).  XI 


(SPHX.SPHY) 
(250).  X9( 250) , 


Y 1(250) , 


Y9( 250) 


3 

4 

5 

6 

7 

8 
9 

10 

1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


C 

C  THIS  SUBROUTINE  CALCULATES  SPROCKET  HOLE  AVERAGES  FOR  EACH  CAMERA 
C  FROM  5  CONTINUOUS  FRAMES  IN  THE  PHOTO  DATA  USING  THE  CHAIR  OR  SLED 
C  DISPLACEMENTS  AND  THE  SPROCKET  HOLE  DISPLACEMENT. 

C  (THESE  ARE  THE  900  ANO  100  RECORDS) 

C 

C  ARGUMENT  DEFINITIONS: 

C  SPHX - RETURNED  SPROCKET  HOLE  AVG.  FOR  X  VALUES 

C  SPHY-RETURNEO  SPROCKET  HOLE  AVG.  FOR  Y  VALUES 

C 

GO  TO  50 
C 
C 

C  SAVE  X  AND  Y  DISPLACEMENTS 
C 

ENTRY  SPH9X  (X) 

DO  10  0-1.250 
10  X9(d)«X(d) 

ICX-ICX+1 

RETURN 

C 

C 

ENTRY  SPH9Y  (X) 

DO  20  0- 1.250 
20  Y9(d)-X(d) 

ICY-ICY+1 

RETURN 

C 

C 

ENTRY  SPH1X  (X) 

DO  30  0-1,250 
30  X1(d)-X(d) 

ICX«ICX*1 

RETURN 

C 

C 

ENTRY  SPH1Y  (X) 

DO  40  0-1,250 
40  Y1(d)-X(d) 

ICY-ICY+1 

RETURN 

C 

C 

50  CONTINUE 

C 

C  CALCULATE  DIFFERENCES  IN  SPROCKET  HOLE  AND  SLED  OR  CHAIR 

C  TARGET  DISPLACEMENTS  FOR  5  CONSECUTIVE  FRAMES. 

C  AND  SUM  THESE  DIFFERENCES  (X  VALUES) 

C 

IBG-0 

IF  (ICX.LT.2)  GO  TO  80 
C 

SUMX-0 

DO  70  0-1.250 

IF  (X1(d). GT. 900. .OR. X9(0).GT. 900. )  GO  TO  60 
SUMX*SUMX+( ABS(X1 (d) )-ABS(X9(d) ) ) 

IF  (IBG.EO.O)  IBG-J 
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60 

61 

IF  ( ( U - I BO ) . GE . 4 )  GO  TO  90 

GO  TO  70 

• 

62 

60 

5UMX-0 

63 

IBG-0 

64 

70 

CONTINUE 

65 

C 

66 

IF  (IBG.GE.250)  IBG-0 

67 

C 

68 

80 

SPHX-0 

A 

69 

GO  TO  100 

70 

90 

CONTINUE 

71 

C 

72 

C 

CALCULATE  SPROCKET  HOLE  AVERAGES  AND  CONVERT  TO  PDS  UNITS 

73 

C 

74 

C 

FROM  CENTIMETERS  (X  VALUE) 

75 

C 

SPHX- ( SUMX/5 . )•( lOOO. /( 1 . 884765- . 0254 ) ) 

76 

_  —  — 

77 

C 

78 

lOO 

IF  (ICY.LT.2)  GO  TO  120 

79 

C 

80 

c 

CALCULATE  DIFFERENCES  IN  SPROCKET  HOLE  AND  SLED  OR  CHAIR 

81 

c 

82 

c 

DISPLACEMENTS  FOR  5  CONSECUTIVE  FRAMES  AND  SUM  THESE  VALUES. 

83 

c 

USE  THE  SAME  5  FRAMES  USED  IN  CALCULATING  THE  X  VALUES  ABOVE 

84 

c 

(V  VALUE) 

... 

85 

c 

86 

IST-4 

87 

IF  (IBG.NE.O)  IST-IBG 

88 

IND-250 

89 

IF  (IBG.NE.O)  I NO- 1ST +4 

90 

c 

91 

SUMY -0 

92 

DO  110  U-IST.INO 

' 

93 

IF  (Y1(J). GT. 900. .OR. V9(U).GT. 900. )  GO  TO  120 

• 

94 

SUMY -SUMY +  ( ABS(  Y  1(0)  ) * ABS(  Y9(  vl)  )  ) 

95 

HO 

CONTINUE 

96 

c 

97 

GO  TO  130 

98 

c 

99 

120 

SPHY-O 

lOO 

GO  TO  140 

tr 

101 

c 

• 

102 

130 

CONTINUE 

*  '■ 

103 

C 

104 

105 

c 

c 

CALCULATE  SPROCKET  HOLE  AVERAGE  AND  CONVERT  TO  PDS  UNITS 

106 

c 

FROM  CENTIMETERS  (Y  VALUE) 

107 

c 

SPHY»( SUMY/5 . )♦( 1000. /( 1 . 884765* . 0254 ) ) 

108 

109 

c 

• 

110 

140 

ICX-0 

111 

ICY-0 

1 12 

c 

113 

114 

1 15 

c 

RETURN 

ENO 

♦PRT.S  . ASC8T9 
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LAMBERT 

«tpf$(o) 

,  ASC8T9( 17) 

1 

SUBROUT INE  ASC8T9( INPUTA . I WO , IB I T , NCHAR , CSTRNG ) 

2 

c 

3 

c 

FUNCTION: 

4 

c 

CONVERT  STRING  OF  8-BIT  ASCII  CHARACTERS  TO 

S 

c 

UNIVAC  COMPATIBLE  9-BIT  ASCII  CHARACTER  STRING 

6 

c 

7 

c 

ARGUMENT  DEFINITIONS: 

8 

c 

INPUTA -INTEGER  INPUT  ARRAV  OF  B-BIT  CHARACTER  DATA 

9 

c 

IWO-STARTING  WORO  IN  INPUTA  OF  8-BIT  ASCII  CHARACTER  STRING 

10 

c 

IBIT-BEGINNING  BIT  IN  IWO  OF  INPUT  8-BIT  ASCII  CHARACTER  STRING 

1 1 

c 

NCHAR -NUMBER  ASCII  CHARACTERS  TO  BE  CONVERTED  TO  UNIVAC  FORMAT 

12 

c 

CSTRNG -RE TURNED  OUTPUT  STRING  OF  UNIVAC  9-BIT  ASCII  CHARACTERS 

13 

c 

14 

c 

PROGRAMMER:  J  LAMBERT  20  MAT  1983 

15 

c 

16 

INTEGER  INPUTA(I) 

17 

c 

18 

c 

LENGTH  OF  CHARACTER  STRING  MUST  BE  INITIALIZED  IN  CALLING  PROGRAM 

19 

c 

20 

CHARACTER* ( * )  CSTRNG , I  TEMP* 1 

21 

DATA  I L 1 , I L2/2*0/ 

22 

IB-IBIT 

23 

IW-IWD 

24 

c 

25 

c 

1ST  BIT  OF  UNIVAC  ASCII  CHARACTER  ALWAYS  O(ZERO) 

26 

c 

27 

BITS( I TEMP , 1 . 1 )*0 

28 

c 

29 

c 

CONVERT  EACH  8-BIT  CHARACTER  TO  9-BITS 

30 

c 

31 

00  20. 1-1. NCHAR 

32 

IF((I8+7).LT. 36 ) THEN 

33 

BITS( I TEMP ,2,8)»BITS( INPUT A ( IW) , IB .8 ) 

34 

ELSE 

35 

IL1-37-IB 

36 

IL2-8-IL1 

37 

BITS( I TEMP  ,2,IL1)-BITS( INPUT A ( IW),IB,IL1) 

38 

BITS(ITEMP.2+IL1.IL2)-BITS(INPUTA(IWM). 1.IL2) 

39 

ENO  IF 

40 

SUBSTRI CSTRNG. I . 1 )=SUBSTR( ITEMP . 1 , 1) 

41 

IB-IB+8 

42 

IF( IB. GT . 36) THEN 

43 

IW-IW+1 

44 

IB-IB-36 

45 

END  IF 

46 

20 

CONTINUE 

47 

RETURN 

48 

ENO 

•PRT.S  . XQTPLT 
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LAMBERT 

1 


■TPF$(0).XQTPLT(11) 

•XOT  MPOO51/PL0T 


[/« 

4 

•  .  POS  5  - 

NO.  INPUT  UNITS  (1.  2  OR  3) 

S 

•  .  POS  lO 

-  TIME  RECORD  ANALYSIS  PLOT  OPTION 

y'' 

6 

9. 

0  -  OUTPUT  TO  LINE  PRINTER 

ys 

7 

9. 

1  -  OUTPUT  TO  COMPRESSED  PLOT  FILE 

M 

8 

9. 

NOT  0  OR  1  •  NO  TIME  RECORD  ANALYSIS 

■ 

9 

9.  POS  15 

-  XY  CONTOUR  PLOT  OPTION 

10 

9. 

0  -  PLOT  XY  CONTOUR  FOR  ALL  MOUNTS/SITES/RUNS 

11 

9. 

NOT  O  •  NO  XY  CONTOUR  PLOTS 

12 

9.  POS  20 

-  FRAME  PLOT  OPTION 

13 

•  . 

0  •  FRAME  PLOTS  FOR  ALL  MOUNTS/SITES/RUNS 

14 

9. 

NOT  0  •  NO  FRAME  PLOTS 

IS 

9.  POS  25 

-  FORMAT  OPTION  FOR  FRAME  PLOTS 

■ 

16 

9. 

0  ■  COL'S  2  ft  4  UP  --  COL'S  1.3*5  DORN 

■ 

17 

9. 

NOT  0  •  ALL  COL'S  DOWN 

18 

9.  POS  35 

-  LOWEST  RUN  NO.  TO  BE  PROCESSED 

19 

•  .  POS  45 

-  HIGHEST  RUN  NO.  TO  BE  PROCESSED 

20 

•XQT  SCC$*1NTLIB . FR80POP 

21 

M040 

22 

DATASYSTEMS 

23 

5-4890 

E 

24 

FRAME  PLOTS 

25 

FRAME  PLOTS 

26 

1  HARDCOPY  1  MICROFICHE 

27 

28 

29 

ORAM* 1 -ENDS 

30 

•OELETE.C  COMPRESS. 

31 

•FREE  COMPRESS. 

32 

•ASG.CP  COMPRESS. .F///1000 

m 

33 

•COPY  17. .COMPRESS. 

34 

•FREE  COMPRESS. 

35 

•.  EXITED  FROM 

. XOTPLT 

•PRT.S 

-PLOT 
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LAMBERT 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


*TPF$!0) 
C 
C 
C 
C 
C 
C 
C 
C 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


PL0T(9) 

FUNCTION: 

THIS  ASCII  FORTRAN  DRIVER  PROGRAM  CONTROLS  THE  GENERATION 
OF  MOUNT  PLOTS  FOR  A  RANGE  OF  RUNS  SPECIFIED  BY  THE  USER. 
TWENTY-FIVE  FRAMES  (OATA  PTS.)  ARE  PLOTTED  ON  EACH  OF  10 
PAGES  FOR  A  GIVEN  MOUNT  AND  SITE  10.  THE  DISSPLA  PLOT 
PACKAGE  IS  USED  FOR  ALL  PLOTTING. 

COMMON/LOAD/XMIN.XMAX.  IPGE.YMIN.YMAX. TITLE.  1ST.  ILT,  IX. XOR. 

* YOR . VARX . VARY . I ERR . STEP , DC ( 250 ) 

COMMON/PHTHOR/RUNIO . VAR ID . I S I T ID . F 1 TIME . JDTD1G. DELTA . IPGVER 
CHARACTERS  RUNID , VARID . JOTDIG , IPGVER . TITLE(9 ) . VARX( 12),VARY( 12) 
CHARACTERS  FRUN. LRUN. RUNS ( 200 ) . FIRST . LAST . ENC*4 . SITE* 1 
DATA  MAX/280/ 

READ  USER  OPTIONS  IN  CARO  IMAGE 

READ(5,6)  IFILE.ITIM.IXY, IMTPLT , I WAY . FRUN. LRUN 
FORMAT ( 515, 9X , 2 ( A6 , 4X ) ) 

WHERE: 

I FILE :  NO.  OF  INPUT  FILES 
ITIM:  TIME  RECORO  ANALYSIS  PLOT  OPTION 
O-OUTPUT  TO  LINE  PRINTER 
1 -OUTPUT  TO  COMPRESSED  PLOT  FILE 
NOT  O  ANO  NOT  1  •  NO  TIME  RECORD  ANALYSIS 
IXY:  XY  CONTOUR  PLOT  OPTION 

O-P LOT  XV  CONTOUR  FOR  EACH  MOUNT/SITE/RUN 
NOT  O  *  NO  XY  CONTOURS 
IMTPLT:  MOUNT  PLOT  OPTION 

O-FRAME  PLOTS  REQUIRED  FOR  ALL  DATA 
NOT  O  *  NO  FRAME  PLOTS 
PLOT  DIRECTION  OPTION 
0- BOTTOM  TO  TOP  COLUMNS  2  AND  4 
NOT  0  -  TOP  TO  BOTTOM  ALL  C0LUW4S 
LOW  ORDER  VALUE  OF  RUN  NO.  RANGE 
HIGH  ORDER  VALUE  OF  RUN  NO.  RANGE 


IWAV : 


FRUN: 

LRUN: 


SET  COMPRESSOR  FOR  PLOTTING 

CALL  COMPRS 
CALL  SET0EV(6,0) 

CALL  GRACE (O. )  * 

CALL  PAGE ( 1 4 . / 1 . 02 . / 1 ■ 02 ) 
CALL  NOBRDR 
CALL  ERASE 

ONE  UNIT  •  10.5/1.02  INCHES 

CALL  UNITS! 10.5/1 .02) 

COMPUTE  RUN  RANGE 

FIRST«SUBSTR( FRUN, 3 , 4 ) 

LAST  *SUBSTR( LRUN ,3,4) 

DECODE (4. FIRST) IF IRST 
DECODE (4. LAST) I LAST 
FORMAT! 14) 
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oO 

C 

61 

C 

COMPUTE  INPUT  UNIT  ASSIGNMENTS 

62 

C 

63 

ILIM-IFILE+8 

64 

IU»9 

65 

1 

CONTINUE 

66 

PRINT  3. IU. ILIM, ITIM, IXY , IMTPLT . 1WAY . FRUN . LRUN 

67 

3 

FORMAT (6I10,2(4X,A6)) 

68 

C 

69 

C 

GENERATE  TITLE  FOR  EACH  PAGE 

70 

C 

71 

TITLE( 2  )» '  CAM' 

72 

TITLE(5 )• ' 

73 

TITLE(6)* '  PAG' 

74 

TITLE ( 7 )  = '  E 

75 

TITLE(9)= ' 

76 

T ITLE( 3  )  = ' 

77 

TITLE (4  )  = ' 

78 

TITLE ( 8  )* ' 

79 

C 

80 

C 

PLOT  DATA  FOR  EACH  RUN 

81 

C 

IRN-I ABS( IL AST- IF  I RST )+ 1 

82 

83 

INC=IFIRST 

84 

DO  500 .1*1,1 RN 

85 

C 

86 

C 

CREATE  RUN  ARRAY  TO  BE  SEARCHED 

87 

C 

88 

SUBSTR(RUNSd).  1.2)-'LX' 

89 

ENCODE ( 17 , ENC ) INC 

90 

17 

FORMAT (14) 

91 

SUBSTR(RUNS(I).3.4)«ENC 

92 

INC” INC* 1 

93 

C 

ATTEMPT  TO  SET  SECTOR  TO  HEADER  FOR  PUN 

94 

C 

95 

C 

CALL  SET5CT( IU.RUNS( I ) . 'PHOATA' , 1 , I SITE ,NREC , IER) 

96 

97 

IF( IER.NE  0)G0  TO  500 

98 

CALL  GETRCD( IU . MAX . RUNID , VARID , I SITE , IDL , DC , I ER ) 
IF ( IER . NE . 0)THEN 

99 

10O 

PRINT  13. IER 

ioi 

13 

FORMAT ( IX , 'STATUS- '.03) 

102 

STOP  'HEADER  ERROR' 

103 

104 

END  IF  A 

IF(VARID.NE.  'PHOATA'  )STOP  'HEADER  ERR®' 

105 

C 

GENERATE  TIME  RECORD  ANALYSIS  ANO  WRITE  ON  FICHE 

106 

c 

107 

c 

108 

ENCODE ( 10,SITE)ISITID 

109 

io 

FORMAT (11) 

110 

SUBSTR( TITLE( 3) ,3, 1)»SITE 

111 

TITLE(  D-RUNIO 

112 

I F ( ITIM.EO.O  .OR.  ITIM.EQ. 1 )THEN 

113 

CALL  TIMER(IU.ITIM) 

114 

IF( IERR.NE .0)00  TO  500 

115 

END  IF 

116 

c 

PLOT  DATA  FOR  EACH  MOUNT  FOR  THIS  RUN 

117 

c 

118 

c 

s 

,  ■-  --  1,-  -  ' 

• 

7." 

■ 

120 

C 

LOAD  DATA 

121 

c 

1 

122 

C 

• 

■ : 

123 

CALL  LDATA(IU.L) 

124 

IF ( IX . NE . 0)THEN 

0 

125 

1ST*  1 

126 

I FM*  1 

127 

ELSE 

s 

128 

GO  TO  400 

m 

12Q 

END  IF 

■ 

130 

c 

-  - 

131 

c 

GENERATE  TITLE  INFORMATION  FOR  THIS  MOUNT 

• 

132 

c 

133 

IF ( L . EQ. 1 )THEN 

134 

TITLE(5)='N£CK  ' 

135 

ELSE  I F ( L . EO . 2 )THEN 

136 

TITLE(5)= 'MOUTH  ' 

■ 

137 

ELSE  IF( L . EO . 3 )THEN 

— 

138 

TITLE(5)*' HEAD  ' 

• 

139 

ELSE  IF( L . EO ■ 4 )THEN 

140 

TITLE(5)=' PELVIS' 

141 

END  IF 

142 

c 

143 

c 

LOOP  FOR  EACH  PAGE  OF  PLOTS  (FRAMES) 

144 

I F( IMTPLT . NE . 0)G0  TO  300 

•*» 

145 

c 

t 

146 

DO  55, IPGE* 1 . 10 

147 

c 

148 

c 

DEFINE  PHYSICAL  ORIGIN  FOR  FIRST  FRAME  ON  PAGE 

149 

c 

150 

XOR-. 16 

151 

Y0R=.81 

1 

152 

c 

i 

153 

c 

DETERMINE  MAX  AND  MIN  VALUES  FOR  THIS  PAGE  OF  PLOTS 

•»  1 

154 

c 

- 

• 

155 

I LT* 1ST +24 

156 

CALL  MIXY 

157 

c 

158 

c 

PRINT  TITLE  AT  TOP  OF  PAGE 

- 

159 

c 

s" 

160 

CALL  PLOT  1 

161 

c 

| 

162 

c 

LOOP  TO  CONTROL  COLUMN  OF  PLOTS 

m 

163 

c 

»E _ _ 

164 

DO  50, K* 1,5 

165 

IC-M0D(K.2) 

166 

c 

167 

c 

PLOT  DATA  5  FRAMES  FOR  EACH  OF  5  COLUMNS 

168 

c 

169 

DO  25, J* 1,5 

| 

170 

CALL  PL0T2OFM) 

' 

171 

c 

W  -  -1 

- 

172 

c 

CHECK  FOR  DIRECTION  OF  PLOTS  (UP  OR  DOWN) 

173 

c 

174 

IF(IWAY.NE.O)THEN 

175 

c 

176 

c 

PLOT  TOP  TO  BOTTOM  ALL  COLUMNS 

177 

c 

m 

178 

IF ( J . NE . 5 ) YOR* YOR- . 2 

m 

179 

ELSE 

1 

i 

K  * 
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180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 


C 

C  PLOT  BOTTOM  TO  TOP  COLUMNS  2  AND  4 
C 

IF(d.NE.5)THEN 

IFUC.NE.O)THEN 

Y0R-Y0R-.2 

ELSE 

Y0R«Y0R+.2 
END  IF 
END  IF 
END  IF 
C 

C  CHECK  FOR  LAST  DATA  POINT 
C 

IFOST.LT. 250)THEN 
IST-IST+1 
IFM*IFM+1 
ELSE 

CALL  ENOPL(O) 

GO  TO  55 
END  IF 

25  CONTINUE 
C 

C  MODIFY  PHYSICAL  ORIGIN  FOR  NEXT  COLUMN 
C 

IF( I WAY .NE . 0)THEN 
XOR-XOR+ . 2 
V0R-.81 
ELSE 

XOR>XOR+ . 2 
END  IF 

50  CONTINUE 

CALL  ENDPL(O) 

55  CONTINUE 
C 

C  PLOT  X  VS  Y  ALL  FRAMES 
C 

300  IF( IXY . EQ. 0)CALL  XYPLOT 
400  CONTINUE 
500  CONTINUE 

IF( IU. LT . ILIM)THEN 
IU-IU+1 
GO  TO  1 
END  IF 
CALL  DONEPL 
STOP 
END 


VPRT.S  . PLOTD 
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6 
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9 
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12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 
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»TPF$!0) .TIMER (33) 

SUBROUTINE  TIMER( IU, ITIM) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

C 

C 


FUNCTION: 

GENERATE 


AND  PLOT  TIME  RECORD  ANALYSIS  FOR  A  GIVEN  RUN 


13 


C 

C 

C 


ARGUMENT  DEFINITIONS: 

IU:  LOGICAL  UNIT  NO.  OF  INPUT  FILE 
ITIM:  OPTION  TO  PLOT  TIME  RECORD  ANALYSIS  WHERE 
0  *  OUTPUT  TO  LINE  PRINTER 
1  -  OUTPUT  TO  COMPRESSED  PLOT  FILE 

COMMON/LOAD/XMIN . XMAX , IPGE , YMIN, YMAX . T I TLE . I ST . I LT , I X , XOR , 

*YOR . VARX . VARY . I ERR . STEP , DC( 250 ) 

COMMON/PHTHDR/RUNID, VARID. I SITID. FITIME. JDTDIG. DELTA, IPGVER 
CHARACTER'S  RUNID , VARID . JDTDIG . IPGVER ,TITLE(9) . VARX( 12) ,VARY( 12) 
CHARACTER'S  SYM* 1 ( 250) . LTINT . LDD*9( 250) . LETIME*9( 250) 

CHARACTER'S  LINDEX*3( 250) . LPDST . LEXPT , LFRN0*7 
CHARACTER* 128  LINE , LHEA0(2) 

DIMENSION  ETIME ( 250) , DD( 250) 

DATA  LPDST. LEXPT. LFRNO/'D  TIME'.'P  TIME*. 'FRAME  '/ 

IERR-1 

CALL  HEIGHT!. 018) 

SET  SECTOR  FOR  TIME  RECORD 
VAR 10*'  TS  ' 

SUB STR( VARID , 1 . 1 ) *SUBSTR( TITLE! 3) . 3 . 1 ) 

CALL  SETSCTCU,  RUNID.  VARID.  1 ,  IS  IT  ,NREC .  IER) 

IF(IER.NE.O) THEN 
PRINT  13. IER 
FORMAT! IX. 'STATUS*'. 03) 

GO  TO  lOO 
END  IF 

READ  TIME  RECORD 
MAX-280 

CALL  GETRCDOU. MAX, RUNIO, VARID. ISIT.IDL.DD. IER) 

IF! ier.ne .0)THEN 
IF(IER.EO. OTHEN 
GO  TO  lOO 
ELSE 

PRINT  13. IER 
STOP  'TIME  ERROR' 

END  IF 
END  IF 

CHECK  FOR  TIME  RECORD  IDENTIFICATION 

IF(SUBSTR! VARID. 2 , 2 ) .NE . 'TS' )THEN 
PRINT  1. RUNID 

FORMAT! IX, 'COULD  NOT  FIND  TIME  RECORD  FOR  RUN  '.A6) 

STOP  'TIME  ERROR' 

ELSE 
IERR*0 
END  IF 

COMPUTE  TIME  INTERVAL  BETWEEN  FRAMES 
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DO  15.K-250. 1,-1 
IF (DO(K ) . LE . 0)THEN 
CONTINUE 
ELSE 

00  21 , d= 1 , 250 
IF (D0( d  )  .  LE . 0)THEN 
CONTINUE 
ELSE 

IF(ODCK)  LT.OD(d) )DD(K)*DD(K)+10. 
TINT»(DO(K)-OD(d) )/(K-d) 

GO  TO  32 
END  IF 
CONTINUE 
GO  TO  100 
END  IF 
CONTINUE 
PRINT  1 6 . RUN I 0 

FORMAT ( IX. 'TIME  RECORO  FOR  RUN  '.A6,'WAS  ALL  ZEROES') 
GO  TO  10O 
ETIME( 1 )-0D( 1 ) 

SYM(1)«'  ' 

COMPUTE  EXPECTEO  TIMES 


DO  35.Ud-2.250 
ETIME(dd)*ETIME(dd-1 )+TINT 
SYM(dd)-'  ' 

IF(OD(dd)  -LT.  . OOOOO 1 ) THEN 
CONTINUE 
ELSE 

DIF*ETIME(dd)-00(dJ) 

I F( ABS(DI F )  .GT.  .  OOOSYMC  dd  )  *  '  • ' 

END  IF 
CONTINUE 

WRITE  HEADER  INFO  FOR  TIME  RECORD  ANALYSIS 

ENCODE( 1 1 , LTINT )TINT 
FORMAT (F6. 5) 

IF(ITIM.EQ. 1)THEN 

CALL  TABLET( 'CENTER' . 'LONG' ) 

ELSE  IF ( ITIM. EO -0)THEN 
PRINT  56 
FORMAT ( ' 1 ' ) 

END  IF 

SUBSTR ( LHE AD (1), 1,4)-' RUN  ' 

SUBSTR(LHEAD( 1 ) .5.6)»RUNI0 
SUBSTR(LHEAO< 0.1 1.11)-'  ' 

SUBSTR(LHEAD( 0.22. 7)- 'CAMERA  ' 

SUBSTRC LHEADC 1 ) . 29 , 3 ) *SUBSTR( TITLE( 3 ) , 1.3) 
SUBSTR (LHE AD( 1) . 32 , 1 O- '  ' 

SUBSTR(LHEAD( 1 ) ,43, 15)*'PREDICTED  RATE  ' 
SUBSTRC  LHEAD( 1 ) ,58,6)-LTINT 
I F  (  ITIM.  EO  .  OTHEN 

SUBSTRC LHE AD ( 1) .64. 1 )- 
CALL  CTLINE( LHEADC 1 ) ) 

CALL  CTLINE( '  $') 

ELSE  IFCITIM.EO. 0)THEN 
PRINT  14 , LHEADC 1 ) 

FORMATC  2 IX . A63/ ) 
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» 


!.'!0 

END  IF 

121 

C 

122 

C 

WRITE  HEADER  INFO.  FOR  INDIVIDUAL  COLUMNS 

123 

C 

124 

CALL  HE  I GHT (.01) 

125 

DO  23, IREP*  1 , 5 

126 

SUBSTR(LHEAD(2).(((IREP-1)»23)+1),7)«LFRN0 

127 

SUBSTR( LHE AD( 2 ) , ( ( ( IREP- 1 ) »23 )+8 ) , 8 )*LPDST 

128 

SUBSTR(LHEAD(2) . ((( IREP- 1 ) *23 )♦ 16 ) . 8)«LEXPT 

129 

23 

CONTINUE 

130 

IF(ITIM.EO. 1)THEN 

131 

SUBSTR( LHEAD( 2 ) , 1 16 , 1)*'$' 

132 

CALL  CTL INE ( LHEAD( 2 ) ) 

133 

CALL  CTLINE( '  $'  ) 

134 

ELSE  IF ( ITIM . EQ . OlTHEN 

135 

PRINT  17 , LHEAD( 2 ) 

136 

17 

FORMAT! IX. A 1 15/) 

137 

END  IF 

138 

C 

139 

C 

PRINT  TIME  RECORD  ANALYSIS 

140 

C 

141 

DO  60 , KK  » 1 , 50 

142 

DO  71,1 ENC= 1 , 5 

143 

INDEX* ( ( < IENC- 1 ) *50)+KK ) 

144 

ENCODE! 8. LINDEX( INDEX)) INDEX 

145 

ENCODE ! 9. LDO! INDEX))  DD(INOEX) 

146 

ENC0DE!9.LETIME( INDEX) )  ETIME! INDEX) 

147 

9 

FORMAT! F9. 5) 

148 

8 

FORMAT! 13) 

149 

SUBSTR!LINE.(!(IENC-1)*23)+1), 1)*'  ' 

150 

SUBSTR(LINE.( ! ! IENC- 1 ) *23 )+2 ) , 3 )»LINDEX( INDEX ) 

151 

SUBSTR! LINE . ! ( ! IENC- 1 ) *23 )+5) , 9)*LDD( INDEX ) 

152 

SUBSTR!LINE. ! ! ! IENC- 1 )*23)+ 14) ,9)=LETIME! INDEX) 

153 

SUBSTR! LINE . ( ! ! IENC- 1 )*23)+23) , 1 )=SYM! INDEX) 

154 

71 

CONTINUE 

155 

IFUTIM.EO.  DTHEN 

156 

SUBSTR! LINE , 116, 1)-'$' 

157 

CALL  CTLINE! LINE ) 

158 

ELSE  IF! ITIM. EO . 0)THEN 

159 

PRINT  18. LINE 

160 

18 

FORMAT! IX, A1 15) 

161 

END  IF 

162 

60 

CONTINUE 

163 

IFUTIM.EO. 0)THEN 

164 

PRINT  78 

165 

78 

FORMAT! ' 1' ) 

166 

END  IF 

167 

ioo 

CONTINUE 

168 

CALL  RESET! 'HEIGHT' ) 

169 

IF!  ITIM.  EO-  DTHEN 

170 

CALL  ENDTAB!0) 

171 

CALL  ENDPL!0) 

172 

END  IF 

173 

RETURN 

174 

END 

•PRT,S 

. XQTSGN 
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LAMBERT* TPFS(O) . LOATA( 13) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


C 

C 

C 

C 

C 

C 

C 

C 

C 


SUBROUTINE  LDATA( IU.MNT ) 
FUNCTION: 

LOAD  PHOTO  DISPLACEMENT 
ACCESS  BY  PLOT  ROUTINES 


DATA  INTO  COMMON  FOR 


OF  INPUT  DATA  FILE 


ARGUMENT  DEF INT IONS : 

IU:  LOGICAL  UNIT  NO. 

MNT:  MOUNT  INOEX 

COMMON/ LOAD/XMIN, XMAX . IPGE . YMIN, YMAX . TITLE , 1ST , ILT , IX.XOR, 

•YOR . VARX . VARY . IERR . STEP , DC ( 250 ) . DX ( 250 . 12 ) . DY<  250. 12 ) 
COMMON/PHTHOR/RUNID, VARID. ISITID.F 1TIME , JDTDIG, DELTA. IPGVER 

CHARACTER*6DRUN?D.VARIO.UDTOIG. IPGVER ,TITLE(9). VARX( 12 ) , VARY( 12 ) 
CHARACTER*2  ONE/ '01 '/ . TGCHK/ '00' / . MTCHK* 1/ '0' / 

CHARACTER*3  TARGET . 0LDRUN*6 
DATA  MAX/280/ 


19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


IX-O 

IY-0 

C 

C  INITIALIZE  CHAIR  DATA  FLAG 

C 

I F ( OLDRUN . NE . RUNID ) I CHA I R-0 
C 

C  SET  SECTOR  TO  TIME  RECORD 

C 

VARID- '  TS  '  , 

SUBSTR( VARID. 1 . 1 )-SUBSTR( TITLE( 3) ,3 , 1) 

CALL  SETSCT( IU. RUNID. VARID. 1 . ISIT.NREC. IER) 
IF(IER.NE.O) THEN 
PRINT  13. IER 

13  FORMAT< IX. 'STATUS-'. 03) 

STOP  ' LDATA ' 

END  IF 


C 

C 

C 

8 


C 

C 

C 


16 

606 


GET  NEXT  DATA  RECORD 

CALL  GETRCD( IU. MAX. RUNID. VARID, ISIT.IDL.DD, IER) 

I F ( I ER . NE . O ) THEN 
IF(IER.EO. OTHEN 
GO  TO  56 
ELSE 

PRINT  13. IER 
GO  TO  72 
END  IF 
END  IF 

CHECK  FOR  CHAIR  DATA  AND  STORE 

IF( I CHAIR . EO. OTHEN 
CONTINUE 
ELSE 

TARGET-SUBSTR(VARID, 3,3) 

IF (TARGET . EO. ' 100' )THEN 
DO  16.d-1.250 
OC(d)«DD(d) 

PRINT  606. VARID, RUNID 

FORMAT ( '  CHAIR  DATA,  VARIABLE ' ,A6, '  LOADED  FOR  RUN  *.A6) 
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23 
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28 

29 
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31 

32 

33 

34 

35 

36 

37 
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41 
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43 
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TPF$<0).MIXY<30) 

SUBROUTINE  MIXY 
C 

C  FUNCTION: 

C  DETERMINE  MINIMUM  AND  MAXIMUM  X  AND  Y  VALUES  FOR 

C  A  PAGE  OF  MOUNT  PLOTS 

C 

COMMON/LOAD/XMIN.XMAX, IPGE ,YMIN, YMAX, TITLE. 1ST. 1LT, IX.XOR, 

* YOR . VARX . VAR Y . I ERR , STEP . OC ( 250 ) . DX ( 250 . 1 2 ) , D Y ( 250 , 1 2 ) 
COMMON/PHTHOR/RUNID . VAR 10 , ISITID. FITIME, UDTD I G . DELT A . I PGVER 
CHARACTER'S  RUNIO . VARID , JDTDIG. IPGVER . TI TLE(9) , VARX( 12).VARY( 12) 
XMIN. YMIN=999. 

XMAX , YMAX'O. 

C 

C  INITIALIZE  COUNTERS 

C 

ISTNEW=IST 

ILTNEW-ILT 

C 

C  CHECK  IF  ALL  PTS  TO  BE  PLOTTED 

C 

I F ( I ERR . EO . 999 ) THEN 
ISTNEW* 1 
ILTNEW250 
END  IF 
C 

C  EXAMINE  OATA  PTS.  FOR  ALL  TARGETS 
C 

DO  20. JK'ISTNEW. ILTNEW 
DO  20, JL* 1 . IX 
DXX* ABS( DX( JK .  JL  )  ) 

DYY“ABS(DY ( JK , JL ) ) 

IF (DXX . LE . 990 .  .AND.  DYY . LE . 990. )THEN 
XMIN* AMI N 1 ( DXX . XMI N ) 

YMIN-AMINI(DYY.YMIN) 

XMAX -AMAX 1 ( DXX . XMAX ) 

YMAX  «  AMAX 1  (  DY  Y  .  YMAX  ) 

END  IF 

20  CONTINUE 
C 

C  EOUATE  X  AND  Y  AXES 
C 

XDIF-XMAX-XMIN 
YDIF'YMAX- YMIN 
IF(XDIF.LE.YDIF )THEN 

XMAX*XMAX+( YOIF-XDIF )/2 . 

XMIN'XMIN- ( YDIF -XOIF )/2 . 

ELSE 

YMAX=YMAX+(XDIF - YDI F )/2 . 

YMIN«YMIN-(XDIF-YDIF)/2. 

END  IF 
C 

C  COMPUTE  STEP  INTERVAL  FOR  PLOT 
C 

IF( IERR .NE . 999 ) THEN 
STEP* ( XMAX -XMIN )/ . 18 
ELSE 

STEP* ( XMAX -XMIN)/ . 8 
END  IF 
RETURN 
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• 


LAMBERT  »TPF$(0) 
1 


PLOTO( 46) 

SUBROUTINE  PLOTO 


2 

C 

3 

C 

FUNCTION: 

4 

c 

CONTROLS  (1)  GENERATION  OF  TITLE  FOR  EACH  PAGE  OF  MOUNT 

5 

c 

PLOTS  ANO  (2)  GENERATION  OF  MOUNT  PLOTS  WITH  DISSPLA 

6 

c 

PLOT  PACKAGE 

7 

c 

8 

COMMON/ LOAD/XMIN . XMAX , I PGE . VM1N . YMAX . TITLE . I ST . I LT . IX , XOR , 

9 

» YOR , VARX , VARY , I ERR , STEP , DC ( 250 ) . OX ( 250 , 12),DY(250, 12) 

10 

COMMON/PHTHOR/RUNID . VARID . ISITID.F 1T1ME, UDTDIG , DELTA . IPGVER 

11 

CHARACTER*6  RUNID , VARID , JDTOIG . I PGVER . TITLE (9 ) , VARX( 12),VARY( 

12 

*.PAG*3 

13 

CHARACTERS  TGT ( 12 ) . TADO. TDROP , TNEW( 12).T0LD( 12 ) . SSCL*8 . SC*7 

14 

CHARACTER*4  FRAM/'  '/.TEMP*2/'  '/ 

IS 

DIMENSION  DAT( 12,2) ,NPT( 12 ) 

16 

c 

17 

c 

PLOT  TITLE  AT  TOP  OF  FRAME 

18 

c 

19 

c 

20 

ENTRY  PLOT  1 

21 

c 

22 

c 

FUNCTION: 

23 

c 

GENERATE  TITLE  INFO.  FOR  PAGE  OF  MOUNT  PLOTS 

24 

c 

25 

CALL  RESET ( 'PHYSOR' ) 

26 

CALL  AREA2D( 1 . , .8) 

27 

CALL  HE I GHT ( .018) 

28 

c 

29 

c 

PRINT  X-SCALE  VALUE  FOR  THIS  PAGE  OF  PLOTS  (METERS) 

30 

c 

31 

SCAL-XMAX-XM1N 

32 

ENCODE ( 1 1 , SC ) SCAL 

33 

11 

FORMAT (F7. 6) 

34 

SUBSTR(SSCL.1.7)-SC 

35 

SUBSTR(SSCL,7, 1 )«'$' 

36 

CALL  MESSAG(SSCL. 100. .39, .925) 

37 

c 

38 

c 

PLOT  TITLE 

39 

c 

40 

ENCODE ( 3 , P AG ) I PGE 

41 

SUBSTR(TITLE(7) ,3.3)-PAG 

42 

3 

FORMAT( 13) 

43 

CALL  HEIGHT (.027) 

44 

CALL  MESSAG( TITLE . 44 , 0 . . .92) 

45 

CALL  ENOGR(O) 

46 

CALL  RESET! 'HEIGHT' ) 

47 

RETURN 

48 

c 

49 

c 

50 

ENTRY  PL0T2(1FM) 

51 

c 

52 

c 

FUNCTION: 

53 

c 

GENERATE  MOUNT  PLOT  FOR  A  FRAME  OF  DATA 

54 

c 

ARGUMENT  DEFINITION: 

55 

c 

IFM:  CURRENT  FRAME  COUNT 

56 

c 

57 

c 

PLOT  DATA  POINTS 

58 

c 

59 

CALL  HEIGHT(.OI) 

dt  J 
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bO 

CALL  PHYSOR(XOR. YOU) 

61 

CALL  AREA2D( . 18. . 18) 

62 

C 

63 

C 

WRITE  BLOCK  NO.  CORRESPONDING  TO  DATA  PT . 

64 

C 

65 

XFM«. 153 

66 

YFM* . 165 

67 

ENCODE (7 , FRAM) IFM 

68 

7 

FORMAT ( 13) 

69 

SUBSTR( FRAM .  4 , 1 )  =  '  $ ' 

70 

CALL  MESSAG( FRAM, 100. XFM , YFM) 

71 

C 

72 

C 

PLOT  SYMBOL  FOR  CHAIR  DATA  (IF  CHAIR  DATA  AVAILABLE) 

73 

c 

74 

IF(DC(IST).GT. 990 ) THEN 

75 

CONTINUE 

76 

ELSE 

77 

XCH= .005 

78 

YCH*. 165 

79 

CALL  MESSAG( 'C' , 1.XCH.YCH) 

80 

ENO  IF 

81 

c 

82 

c 

SCALE  DATA  TO  PRESENT  PLOTTING  AREA 

83 

c 

84 

CALL  GRAF ( XMIN. STEP . XMAX , YMIN, STEP . YMAX ) 

85 

c 

86 

c 

DRAW  BLOCK  AROUND  SUBPLOT 

87 

c 

88 

CALL  RLVEC( XMIN. YMIN. XMIN. YMAX, 0000) 

89 

CALL  RLVEC ( XMAX , YMAX . XMIN . YMAX ,0000) 

90 

CALL  RLVEC ( XMAX . YMAX , XMAX . YMIN . OOOO ) 

91 

CALL  RLVEC (XMIN, YMIN, XMAX. YMIN. 0000) 

92 

NTHIS=0 

93 

DO  50,  d* 1 , IX 

94 

NPT(J)-0 

95 

I F  (  DX  (  I  ST .  d  )  .  GT .  990 . OR . OY ( I ST . d ) . GT .  990  )  THEN 

96 

CONTINUE 

97 

ELSE 

98 

NTHIS=NTHIS+1 

99 

c 

lOO 

c 

LOAD  DATA  FOR  ALL  TARGETS  IN  THIS  FRAME 

101 

c 

102 

DAT ( NTHIS . 1 ) *ABS(OX( 1ST , d) ) 

103 

DAT(NTHIS.2)*ABS(0Y(IST,d)) 

104 

c 

105 

c 

IDENTIFY  MANUAL  ENTRIES 

106 

c 

107 

IF (DX( 1ST , d ) . LT .0 . )NPT(NTHIS )* 1 

108 

TGT ( NTHIS ) * SUBSTR ( VARX(d)  ,4,2) 

109 

I F ( SUBSTR ( VARX ( d ) , 4 , 2 ) . NE . SUBSTR ( VAR Y ( d ) , 4 , 2 ) ) THEN 

110 

PRINT  52 

111 

52 

FORMAT ( '  X  +  Y  DISP.  TARGETS  00  NOT  AGREE  ') 

1 12 

STOP  'PLOTD' 

113 

ENO  IF 

1 14 

END  IF 

115 

50 

CONTINUE 

116 

IF(NTHIS.NE.O)THEN 

117 

DO  60, d* 1 .NTHIS 

118 

X«XP05N(DAT ( d , 1 ) . DAT ( d , 2 ) ) 

119 

Y«YPOSN(DAT(d. 1).DAT(d,2)) 

120 

IF(TGT(d) .  LE  .  '09'  )THEN 

121 

XT  -X 

122 

YT«Y 

123 

XM*X- .004 

124 

YM«Y- .009 

125 

SUBSTR(TEMP, 1 , 1 )*SUBSTR( TGT(d) ,2.1) 

126 

IT'  1 

127 

ELSE 

128 

XT*X- . 006 

129 

YT»Y 

130 

XM-X- . 007 

131 

YM-Y - . 009 

132 

TEMP-TGT(d) 

133 

IT-2 

134 

END  IF 

135 

C 

136 

C 

PLOT  TARGET 

137 

C 

138 

CALL  M£SSAG(TEMP, IT .XT. YT) 

139 

C 

140 

C 

PLOT  CIRCLE  AROUNO  TARGET  IF  MANUAL  ENTRY 

141 

C 

142 

IF(NPT(d).NE.O)THEN 

143 

CALL  HEIGHT ( . 029) 

144 

CALL  INTNO(O.XM.YM) 

145 

CALL  HEIGHT ( .01 ) 

146 

END  IF 

147 

60 

CONTINUE 

148 

INS»NTHIS-1 

149 

C 

150 

C 

CONNECT  TARGETS  WITH  LINES 

151 

C 

152 

DO  68.d-1.INS 

153 

Jf»d+1 

154 

DO  68. d2=d1. NTHIS 

155 

CALL  RLVEC(OAT(d. 1 ) . DAT< d , 2 ) ,DAT( d2 , 1).DAT(d2,2) 

156 

68 

CONTINUE 

157 

END  IF 

158 

TADD- '  ' 

159 

TDROP- '  ' 

160 

C 

161 

C 

DETERMINE  WHICH  TARGETS  HAVE  BEEN  ADDED  OR  DELETED 

162 

C 

163 

IF( IST.NE . 1)THEN 

164 

IF(NTHIS.NE.O.OR.NLAST.NE.O)THEN 

165 

IF(NTHIS.NE.O)THEN 

166 

DO  1 12 , d* 1 .NTHIS 

167 

112 

TNEW( d)-TGT(d) 

168 

I F ( NL AST . NE . 0 ) THEN 

169 

DO  1 13 , d- 1 .NLAST 

170 

DO  1 13. K-1 .NTHIS 

171 

IBOTH-O 

172 

I F ( TOLD( d ) . EO . TNEW(K) ) I BOTH* 1 

173 

IF(IBOTH.EO. DTHEN 

174 

TOLD(d)*'  ' 

175 

TNEW(K)*'  ' 

176 

END  IF 

177 

113 

CONTINUE 

178 

DO  114, d»1, NLAST 

179 

114 

IF(TOLO(d) .NE. '  ' )TOROP-TOLD(d) 

L-4 


L  J 


L 
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l 


180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 


ENO  IF 

00  1 15. K»1 .NTHIS 

115  IF(TNEW(K) .NE. '  ' )TADD*TNEW(K) 

ELSE  IF(NTHIS.NE.O)THEN 

IF (NLA ST . EQ. 0)THEN 
DO  214,d=1.NLAST 

214  I F( T0LD( J ) . NE . '  ' )TDROP«TOLD(d) 

ENO  IF 
ENO  IF 

ELSE  IF(NTHIS . EO .0  .AND.  NLAST.EQ.0)THEN 
GO  TO  I  If 
END  IF 

ELSE  IF ( 1ST  .EO.  1 )THEN 
I F(NTHI S . EQ . 0)G0  TO  111 
ENO  IF 
C 

C  STORE  TARGET  NOS.  FOR  THIS  FRAME 
C 

00  116 .0-1. NTHIS 

116  TOLD(d)*TGT(d) 

NLAST*NTHIS 

C 

C  PLOT  SYMBOL  IF  NEW  TARGET  WAS  VISIBLE  IN  THIS  FRAME 
C 

XNEW* . 005 
YNEW-.007 

I F (TAOO . LE . '09' )THEN 

SUBSTR( TEMP , 1. 1 )=SUBSTR(TADD . 2 . 1) 

IA»1 

ELSE 

TEMP-TADO 
IA-2 
ENO  IF 

CALL  MESSAG(TEMP.IA.XNEW.YNEW) 

C 

C  PLOT  SYMBOL  IF  A  TARGET  WAS  OROPPED  IN  THIS  FRAME 
C 

XDROP*. 159 
YDR0P-.0O7 

I F ( SUBSTR ( TOROP , 1 . 1).EQ. 'O' )SUBSTR(TDROP . 1,1)='  ' 
CALL  MESSAG( TOROP . 2 , XOROP , YOROP ) 

111  CONTINUE 
C 

C  ENO  SUBPLOT 

C 

CALL  RESET( 'HEIGHT' ) 

CALL  ENOGR(O) 

RETURN 

END 
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2 

3 

4 

5 

6 

7 
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9 
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11 

12 

13 
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15 

16 

17 
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19 

20 
21 
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23 

24 

25 

26 

27 

28 

29 
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31 

32 

33 

34 

35 

36 

37 

38 

39 

40 
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1  f PF$(0) . XYPLOT ( 22 ) 

SUBROUTINE  XYPLOT 


C 

C 

C 

C 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 


FUNCTION: 

PLOT  X  AND  Y  DISPLACEMENT  VALUES  FOR  GIVEN  RUN  AND  ANATOMICAL 
MOUNT;  EACH  PHOTO  TARGET  PLOTTED  ON  SAME  GRAPH 

DIMENSION  0XARY(250) ,0YARY(25O) 

COMMON/LOAD/XMIN, XMAX , IPGE , YMIN, YMAX . TI TLE , I ST . I LT . IX . XOR . 

•YOR, VARX, VARY, IERR, STEP, DC( 250) . DX( 250. 12).0Y(250, 12) 
COMMON/PHTHDR/RUNIO . VARID . I SITID, F 1TIME , JDTDIG , DELTA . IPGVER 
CHARACTER'S  RUNID , VARID , JDTDIG. IPGVER . TI TLE ( 9) , VARX( 12),VARY(12) 
CHARACTER'S  XYTIT<9) ,TAR*2 ,TIME*8 .DATE *8 
DATA  NPTS, ITAR/2*0/ 

LOAD  TITLE  INFO. 

DO  12,1*1,9 
XYTIT( I)*TITLE(I) 

12  CONTINUE 

XYTIT(5)* '  ' 

XYTITI6)-'  ' 

XYTIT( 7  )  * '  ' 

XYTIT(8)*' 

XYTIT(9)*'  ' 

CALL  ADATE ( DAT E , T IME ) 

XYTIT(4)*TITLE(5) 

SUBSTR(XYT1T(5),4.2)=SUBSTR(DATE. 1.2) 

SUBSTR(XYTIT(6). 1 . 2 ) -SUBSTR (DATE , 3 , 2 ) 

SUBSTR(XYTIT(6) , 4 , 2 )*SUBSTR(0AT£ .5.2) 

SUPPRESS  LISTING  FOR  OUT  OF  RANGE  POINTS 

CALL  NOCHEK 

CALL  RESET ( ' PHYSOR ' ) 

LABEL  AXES 

CALL  AREA2D( .8. .8) 

CALL  XNAME ( '  FILM  PLANE  X  DISPLACEMENT  (METERS) '. 35) 

CALL  YNAME ( '  FILM  PLANE  Y  DISPLACEMENT  (METERS )'. 35 ) 

PRINT  HEADING 

CALL  HEIGHT ( . 027 ) 

CALL  MESSAG(XVTIT ,50,0. . . 92 ) 

CALL  HEIGHT( .016) 

COMPUTE  MIN  AND  MAX  VALUES  FOR  X  AND  Y 

I ERR *999 
CALL  MIXY 
IERR-0 

SCALE  DATA  FOR  PLOTS 

CALL  GRAF ( XMIN, STEP , XMAX , YMIN , STEP . YMAX ) 

PLOT  X-Y  CONTOUR  FOR  EACH  TARGET 


•  J 


JL-J 
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60 

CALL  HEIGHT (.022) 

61 

DO  10.I-1.IX 

62 

LINE-MODU  .5) 

63 

IF (LINE . EQ.O)THEN 

64 

CALL  DOT 

65 

ELSE  IF(LINE.EQ. 1)THEN 

66 

CALL  RESET('DOT') 

67 

ELSE  IF(LINE . EQ. 2 )THEN 

68 

CALL  CHNDOT 

69 

ELSE  IF(LINE.EQ.3)THEN 

70 

CALL  OASH 

71 

ELSE  IF ( LINE . EQ. 4)THEN 

72 

CALL  CHNDSH 

73 

END  IF 

74 

C 

75 

C 

LOAD  OATA  INTO  PLOT  ARRAY  AND  CHANGE  NEGATIVE 

76 

C 

77 

00  30. d* 1.250 

78 

OXARY (  d  ) -DX (  d ,  I ) 

79 

OYARY(d)«DV( d.I ) 

80 

IF (OXARY ( d ) . LT . 0. )THEN 

81 

DXARY (d ) » -  1 . -DXARY (d ) 

82 

END  IF 

83 

IF(OYARY ( d) . LT . 0. )THEN 

84 

OVARY ( d  )  *  -  1 . *OYARY ( d ) 

85 

END  IF 

86 

30 

CONTINUE 

87 

C 

88 

C 

PLOT  DATA  (ONLY  NON-ZERO  PTS) 

89 

C 

90 

IPOS-1 

91 

ITAR-0 

92 

5 

DO  35. d- I PCS. 250 

93 

IF(DXARY(d) . GT . 990. . OR.DYARY ( d) . GT , 990. )THEN 

94 

CONTINUE 

95 

ELSE 

96 

NPTS-O 

97 

IPLOT-d 

98 

GO  TO  15 

99 

END  IF 

10O 

35 

CONTINUE 

101 

GO  TO  10 

102 

15 

00  40, K- I PLOT , 250 

103 

IF(DXARY(K).GT. 990 . . OR.DYARY(K) . GT.990. )THEN 

104 

IF( ITAR.E0-O)THEN 

105 

SX-XPOSN( DXARY ( I PLOT ) ,OYARY( IPLOT ) ) 

106 

SY- YPOSN( DXARY ( IPLOT) ,OYARY( IPLOT) ) 

107 

TAR-SUBSTR( VARX( I ) , 4 , 2) 

108 

IF(SUBSTR(TAR, 1 , 1 ) . EO. '0' )SUBSTR(TAR . 1 , 1 ) • 

109 

CALL  MESSAG(TAR,2,SX,SY) 

110 

ITAR-1 

111 

END  IF 

112 

CALL  CURVE (DXARY( IPLOT ).DYARY( IPLOT ),NPTS.O) 

113 

IPOS-K+ 1 

114 

IF(K. LT . 250) GO  TO  5 

115 

ELSE 

116 

NPTS-NPTS+1 

1 17 

IF(K . EQ. 250)THEN 

118 

IF( ITAR . EO.O)THEN 

119 

SX-XPOSN( OXARY ( I PLOT ) . OY ARY (IPLOT)) 

SIGNS 


*  t 


1 


71 


120 

S V  « YP0SN( DXARY ( JPLOT ) . OVARY ( IPLOT ) ) 

121 

TAR«SUBSTR(VARX( I ) .4,2) 

122 

IF(SUBSTR(TAR, 1 . 1 ) .£0. 'O' )SUBSTR(TAR. 

123 

CALL  MESSA6(TAR.2.SX.SY) 

124 

ITAR* 1 

125 

END  IF 

126 

CALL  CURVE ( DXARY ( IPLOT) ,DYARY( IPLOT) .NPTS 

127 

END  IF 

128 

END  IF 

129 

40 

CONT INUE 

130 

10 

CONTINUE 

131 

CALL  RESET ( 'HEIGHT' ) 

132 

CALL  RESET ( 'XNAME' ) 

133 

CALL  RESE  T ( ' YNAME ' ) 

134 

CALL  RESET ( 'DOT ' ) 

135 

CALL  ENOPL(O) 

136 

RETURN 

137 

END 

•PUT . S  .TIMER 


LAMBERT «TPF$(0) . XQTSGN( 3) 

1  PXQT  MP005 1/SIGN 

2  O 

3  P.  CONTROL  STATEMENT  IMMEDIATELY  AFTER  PXQT  STATEMENT 

4  P.  POS  5  -  NO.  INPUT  UNITS  (1.  2  OR  3) 

5  P.  POS  10  -  LOWEST  RUN  NO.  TO  BE  PROCESSED 

6  P.  POS  20  -  HIGHEST  RUN  NO.  TO  BE  PROCESSED 

7  PTEST  TNE/0/S6 

8  PJUMP  OUT 

9  PASG.TF  PCSXXXXXXXX. .U9S.SAVE05  .  PCSXXXXXXXX 

10  PCOPY.GM  9. .PCSXXXXXXXX. 

11  PCOPY.GM  9. .PCSXXXXXXXX. 

12  PCOPY.GM  9. .PCSXXXXXXXX. 

13  PTEST2: 

14  PTEST  TG/1/S6 

15  PJUMP  TEST3 

16  PCOPY.GM  10. .PCSXXXXXXXX. 

17  PCOPY.GM  10. .PCSXXXXXXXX. 

18  PCOPY.GM  10. .PCSXXXXXXXX. 

19  PTEST3 : 

20  PTEST  TG/2/S6 

2 1  PJUMP  OUT 

22  PCOPY.GM  11. .PCSXXXXXXXX. 

23  PCOPY.GM  1 1 . .PCSXXXXXXXX. 

24  PCOPY.GM  1 1 . .PCSXXXXXXXX. 

25  POUT : 

26  PFREE.I  PCSXXXXXXXX. 

27  P.  EXITED  FROM  . XQTSGN 


PPRT.S  . SGFLIP 


LAMBERT *  TPF$(0) .  SGFLIP(6 1  ) 

1  C 

2  C  FUNCTION: 

3  C  THIS  ASCII  FORTRAN  DRIVER  PROGRAM  PERFORMS 

4  C  (1)  SIGN  CHANGE  OF  NEGATIVE  DATA  PTS.  (MANUAL  ENTRIES) 

5  C  TO  POSITIVE  AND  (2)  WRITES  REFORMATTED  DATA  RECORD  TO  MASS  STORAGE 

6  C 

7  DIMENSION  DC(250) 

8  COMMON/PHTHDR/RUN ID , VAR I 0 . I S I T I D . F 1 T I ME , JDTD 1 G . DE  LT A , I PGVER 

9  CHARACTER’S  RUNID, VARID , JDTDIG, I PGVER , FRUN, LRUN, RUNS (200) 


10 

*,  FIRST, LAST. ENC*4, SITE*  1 

11 

C 

12 

c 

READ  USER  PARAMETERS  IN  CARD  IMAGE 

13 

c 

14 

RE AD (5, 6) I  FILE, FRUN, LRUN 

IS 

6 

FORMAT ( 15 , 9X . 2( A6 , 4X ) ) 

16 

DATA  MAX/320/ 

17 

c 

18 

C 

COMPUTE  RUN  RANGE 

19 

C 

20 

FIRST*SUBSTR( FRUN, 3 , 4 ) 

21 

LAST  *SU8STR( LRUN, 3,4) 

22 

DECODE (4. FIRST) IF IRST 

23 

DECODE ( 4 , LAST ) I  LAST 

24 

4 

FORMAT ( 14 ) 

25 

C 

26 

c 

COMPUTE  INPUT  UNIT  ASSIGNMENTS 

27 

c 

28 

ILIM*IFILE+8 

29 

IU=9 

30 

1 

CONTINUE 

31 

PRINT  3. IU, ILIM, FRUN, LRUN 

32 

3 

FORMAT (2110, 2( 4X , A6 ) ) 

33 

NIU*~ 1*IU 

34 

c 

35 

c 

BEGIN  RUN  LOOP 

36 

c 

37 

IRN*I ABS( I LAST- 1 FIRST  )* 1 

38 

INC-IFIRST 

39 

00  100, 1*1, IRN 

40 

SUBSTR(RUNSd).  1,2)*'LX' 

41 

ENCODE ( 1 7.ENC)INC 

42 

17 

FORMAT ( 14) 

43 

SUBSTR(RUNS(I),3,4)«ENC 

44 

INC* INC* 1 

45 

c 

46 

C 

ATTEMPT  TO  SET  SECTOR  TO  HEADER  FOR  RUN 

47 

C 

48 

CALL  SETSCTdU.RUNS(I), 'PHOATA' . 1 , ISITE .NREC , IER) 

49 

IF(IER.NE.O)GO  TO  100 

50 

CALL  GETRCD(IU,MAX,RUNIO,VARID. ISITE. IDL, DC, IER) 

51 

IF( IER.NE.O)THEN 

52 

PRINT  13, IER 

53 

13 

FORMAT( IX, 'STATUS*' .03) 

54 

STOP  'HEADER  ERROR' 

55 

END  IF 

56 

C 

57 

C 

SET  SECTOR  FOR  TIME  RECORO 

58 

C 

59 

VARIO* '  TS 
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ENCODE (11,SITE)ISIT1D 
FORMA T( II) 

SUBSTR! VARIO , 1 . 1 )*SITE 

CALL  SETSCT ( IU , RUN1D , VARID , 1 . 1  SI T . NREC . I ER ) 

IF(IER.NE.O)GO  TO  110 

READ  NEXT  RECORD 

DO  10. K* 1 , lOO 
CALL  SECTOR( IU„ ISEC) 

CALL  GETRCO! IU , MAX . RUNID . VARIO . ISIT . IDL . DC .  I ER ) 

IF( IER.NE . 0)THEN 
IF( IER.EO.  1 ) THEN 
GO  TO  lOO 
ELSE 

GO  TO  110 
END  IF 
END  IF 
INEG=0 

CHANGE  SIGNS  OF  NEGATIVE  NOS. 

DO  8 ,M= 1 , 250 
IF(DC(M).LT.O)THEN 
INEG= 1 

DC(M)*-DC(M) 

Ed®  IF 
CONTINUE 

WRITE  RECORD  BACK  TO  MASS  STORAGE 

IF( INEG. NE . 0)THEN 

CALL  SET ADR ( IU , ISEC ) 

CALL  PUTRCO(NIU, RUNID. VARID. ISIT, IDL. DC. IER) 
IF(IER.NE.O)GO  TO  110 
END  IF 
CONTINUE 

ERROR  ROUTINE 

PRINT  111, RUNID, IER 

FORMAT!'  ERROR  IN  SIGN  FLIP  ROUTINE  FOR  RUN  ',A6/'  STATUS 
CONTINUE 

I F ( IU . LT . I LIM )THEN 
IU*IU+1 
GO  TO  1 
END  IF 


SUBROUTINE  SETSCT 


a.  Function:  Subroutine  SETSCT  Is  an  ASCII  FORTRAN  program  designed  to 
control  the  po’sTt  1 on 1 ng  of  sector  address  for  either  a  photo  header  record  or 
a  photo  time  record.  The  calling  program  must  provide  the  logical  unit  number 
of  the  photo  data  file  and  the  ASCII  variable  RUNID;  these  arguments  are  used 
to  locate  the  desired  run  in  the  random  access  directory.  In  addition,  the 
user  provides  the  variable  name  (see  Calling  Sequence)  and  the  desired  occur¬ 
rence  of  the  time  record  as  arguments.  There  may  be  more  than  one  POS  (Photo 
Digitizing  System)  physical  record/target  in  a  run,  hence  it  is  necessary  to 
designate  which  one  Is  desired  (even  If  there  is  only  one  PDS  physical  record/ 
target) .* 

Entry  DATSCT  of  subroutine  DIRBFR  Is  called  in  order  to  retrieve  the  sector 
address  of  the  desired  record.  This  address  Is  used  in  a  subsequent  cal)  to 
Entry  SETADR  of  subroutine  FORTIO  in  order  to  correctly  position  the  data 
file. 

This  subroutine  returns  the  site  ID  of  the  data  file  as  argument  ISITE. 
Additionally,  the  number  of  PDS  physical  records  required  per  target  for  the 
desired  run  is  returned  In  argument  NREC. 

b.  Subroutines  (see  reference  1): 

DIRBFR/Entry  DATSCT,  FORTIO/Entry  SETADR 

c.  Calling  Sequence: 

CALL  SETSCT  (IUNIT,  RUNID,  VARID,  NPHYS,  ISITE,  NREC,  IER) 


*  In  the  event  more  than  250  frames  of  photo  data  were  digitized  for  a 
target,  an  extra  PDS  physical  record  Is  required  for  each  additional  250 
frames. 
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SUBROUTINE  SETSCT  (Continued) 


where: 


ARGUMENT 

DESCRIPTION 

TYPE 

IUNIT 

Logical  unit  number  of  photo  directory 
and  data  file 

Integer 

RUNIO 

Run  number  for  requested  sector 
positioning 

6  char 

VARID* 

Variable  name  for  requested  sector 
positioning 

6  char 

NPHYS 

Desired  occurrence  of  time  record 
(Nphys=l  unless  more  than  250  frames 
were  digitized) 

Integer 

I  SITE 

Returned  site  ID  of  data  file 

Integer 

NREC 

Returned  number  PDS  physical  records 
required  for  each  photo  target 

Integer 

IER 

Returned  error  status  code  (See 

Appendix  4e  for  code  definitions) 

Octal 

*  This  argument  Identifies  whether  the  sector  Is  to  be  set  for  a 
header  record  (VARID^PHOATA1 )  or  a  time  record  (VARID='1TS  1 
for  site  1). 


A 


A 
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SUBROUTINE  GETRCD 


a.  Function:  This  subroutine  Is  an  ASCII  FORTRAN  program  designed  to 
control  the  execution  of  sequential  binary  reads  of  FORTRAN  V  photo  (1)  header 
records,  (2)  time  records  and  (3)  displacement  records.  GETRCD  utilizes  Entry 
READIO  of  subroutine  FORTIO  to  retrieve  the  data  record  from  mass  storage  (the 
data  record  is  stored  In  named  common  IOBUFF).  It  then  extracts  the  run  num¬ 
ber  (RUNID)  and  variable  name  (VARIO)  from  the  data  record  and  converts  each 
from  FIELDATA  to  ASCII  with  system  function  subprogram  FD2ASC;  these  converted 
character  variables  are  returned  to  the  calling  program.  The  site  ID  ( I  SITE) 
of  the  data  record  Is  also  returned  as  an  argument. 

If  the  retrieved  record  Is  either  a  time  record  or  a  displacement  record,  the 
word  length  of  the  data  array  (IDL)  and  the  data  array  itself  (DATARY)  are 
returned  to  the  calling  program.  If  the  retrieved  record  is  a  header  these 
arguments  have  no  meaning.  However,  when  a  header  record  is  read  GETRCD  calls 
subroutine  DECHED  In  order  to  (1)  convert  all  alphanumeric  data  In  the  header 
to  ASCII  and  (2)  store  all  header  information  in  named  common  PHTHDR. 

The  user  specifies  MAXLEN  as  the  maximum  word  length  of  the  data  record  being 
read  from  mass  storage.  GETRCD  utilizes  function  subprogram  EXPAND  to  Insure 
that  the  user  has  specified  a  maximum  length  large  enough  to  accommodate  those 
sentinels  and  checksums  which  will  be  read  from  mass  storage  as  part  of  the 
Input  record.  If  MAXLEN  Is  too  small,  an  appropriate  error  status  code  is 
returned. 

b.  Subroutines  (see  reference  1): 

FORTIO/Entry  READIO,  FD2ASC,  ASC2FD ,  DECHED,  EXPAND 


c.  Calling  Sequence: 

Call  GETRCD  (IUNIT,  MAXLEN,  RUNID,  VARID,  ISITE,  IDL,  DATARY,  IER) 
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SUBROUTINE  GETRCD  (Continued) 


where : 


ARGUMENT 

DESCRIPTION 

TYPE 

IUNIT 

Logical  unit  number  of  photo  data  and 
directory  file 

Integer 

MAXLEN 

Maximum  word  length  of  input  record 

Integer 

RUNID 

Returned  run  number  of  input  record 

6  char 

VARID* 

Returned  variable  name  of  input  record 

6  char 

I  SITE 

Returned  site  ID 

Integer 

IDL 

Returned  word  length  of  data  array 

Integer 

DATARY 

Returned  data  array 

Real 

IER 

Returned  error  status  code  (see 

Appendix  4e  for  code  definitions) 

Octal 

d.  Common  Block  Usage: 

COMMON/ IOBUFF/ I 0BUFF{ 4100) 

where: 

VARIABLE 

DESCRIPTION 

TYPE 

I0BUFF(4100) 

Input  array  containing  FORTRAN  (FIELDATA) 
photo  data  record  read  from  mass  storage 
(see  Figures  2c-2e  for  data  record  formats 

Integer 

The  contents  of  this  array  are  passed  as  an  argument  to  the  header  decoding 

subroutine  (DECHED).  Alternatively,  if  the  photo  input  record  is  either  a 
time  record  or  a  displacement  record,  the  data  array  is  extracted  from  IOBUFF 
for  return  to  the  calling  program. 


VARIO='PHDATA'  for  a  photo  header  record 

VARID='ITS  '2TS  '  or  ' 3TS  '  for  a  photo  time  record  (first 
character  represents  site  ID) 

See  Figure  4  for  possible  values  of  YARID  for  a  photo  displacement  record 
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SUBROUTINE  PUTRCD 


a.  Function;  Subroutine  PUTRCD  Is  an  ASCII  FORTRAN  program  designed  to 
control  the  execution  of  sequential  binary  writes  of  FORTRAN  (FIELDATA)  photo 
(1)  header  records,  (2)  time  records,  and  (3)  displacement  records.  This  rou¬ 
tine  provides  the  user  with  complete  directory  maintenance.  The  user  may  uti¬ 
lize  PUTRCD  to  (1)  create  a  photo  data  output  file  with  accompanying  directory 
or  (2)  perform  sequential  output  of  photo  records  without  directory  mainte¬ 
nance.  If  the  user  desires  to  perform  sequential  writes  without  building  a 
corresponding  directory,  (s)he  must  transmit  the  negative  of  the  logical  unit 
number  (of  the  photo  data  file)  as  the  value  of  the  argument  IUNIT  in  the  call 
statement. 

In  the  event  a  photo  header  record  is  to  be  written,  it  is  the  user's  respon¬ 
sibility  to  provide  the  ASCII  version  of  the  header  In  named  common  PHTHDR 
(see  Figure  2c  for  format).  However,  If  the  user  utilizes  input  subroutine 
GETRCD  to  read  an  existing  header  record.  Its  ASCII  version  is  automatically 
stored  In  common.  This  ASCII  version  of  the  header  must  be  converted  to 
FIELDATA;  this  Is  accomplished  by  a  call  to  subroutine  ENCHED. 

The  user  must  provide  PURTCD  with  the  (1)  logical  unit  number  of  the  photo 
data  output  file  (IUNIT),  (2)  run  number  of  the  output  record  (RUNID),  (3) 
variable  name  of  the  output  record  (VARID*),  (4)  site  ID  of  the  output  record 
( I SITE) ,  (5)  word  length  of  the  output  data  array  (IDL)  and  (6)  output  data 
array  (DATARY).  Values  for  items  (5)  and  (6)  must  be  assigned  for  all  output 
records  except  header  records.  Items  (1)  through  (4)  are  passed  by  PUTRCD  to 
subroutine  DIRCTY  where  the  necessary  directory  entries  are  created. 

Once  the  directory  has  been  updated,  PUTRCD  uses  the  word  length  of  the  output 
array  (IDL)  and  the  output  array  Itself  (DATARY)  to  construct  the  output 
record.  Entry  VIRITIO  of  subroutine  FORTIO  is  then  called  to  perform  the 
sequential  write.  The  output  record  Is  passed  down  In  an  argument  array 
( IOBUFF) ,  and  Is  also  stored  In  named  common  IOBUFF. 

b.  Subroutines  (see  reference  1): 

ASC2FD ,  ENCHED,  FORTIO/Entry  WRITIO 

c.  Calling  Sequence: 

Call  PUTRCD  (IUNIT,  RUNID,  VARID,  ISITE,  IDL,  DATARY,  IER) 


VARID* 'PHDATA'  for  a  photo  header  record 

VARID*' ITS  ',  '2TS  '  or  '3TS  '  for  a  photo  time  record  (first 
character  represents  site  ID) 

See  Figure  4  for  possible  values  of  VARID  for  a  photo  displacement  record 
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SUBROUTINE  PUTRCD  (Continued) 


where: 


ARGUMENT 

DESCRIPTION 

TYPE 

IUNIT 

Logical  unit  number  of  photo  data 
and  directory  file 

Integer 

RUNID 

Run  number  of  output  data  record 

6  char 

VARID 

Variable  name  of  output  data  record 

6  char 

I  SITE 

Site  ID  of  output  record 

Integer 

IDL* 

Word  length  of  output  data  array 

Integer 

DATARY* 

Output  data  array 

Real 

IER 

Returned  error  status  code  (see 

Appendix  4e  for  code  definitions) 

Octal 

d.  Common  Block  Usage: 

(1)  COMMON/PHTHDR/RUNID,  VARID,  ISITID,  FITIME,  JDTDIG,  DELTA,  IPGVER 

(See  Figure  2c  for  format  of  named  common  PHTHDR) 

(2)  COMMON/I OBUFF/IOBUFF( 4100) 

where: 

VARIABLE  DESCRIPTION  TYPE 

I0BUFF(4100)  Output  array  containing  FORTRAN  (FIELDATA)  Integer 
photo  (1)  header  record,  (2)  time  record 
or  (3)  displacement  record  to  be  written 
to  core  (see  Figures  2c-2e  for  data  record 
formats) 

*  This  argument  has  no  meaning  for  header  record. 
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SUBROUTINE  EOFPUT 


a.  Function:  This  ASCII  FORTRAN  program  Is  designed  to  control  the 
writing  of  an  "end-of-flle"  to  the  end  of  a  data  file.  Initially  EOFPUT  calls 
subroutine  DIRCTY  in  order  to  have  the  output  file  positioned  for  the  write. 
Once  this  Is  done,  the  "end-of-flle"  Is  written  by  virtue  of  a  call  to  entry 
WRITIO  of  the  subroutine  FORTIO.  Control  is  then  returned  to  the  calling  pro¬ 
gram.  In  the  event  the  user  wishes  to  write  an  “end-of-flle"  without  direc¬ 
tory  maintenance,  {s)he  must  pass  down  the  negative  of  the  logical  unit  number 
(of  the  photo  data  file)  as  the  value  of  the  argument  IUNIT. 

b.  Subroutines  (see  reference  1): 

DIRCTY,  FORTIO/Entry  WRITIO 

c.  Calling  Sequence: 

Call  EOFPUT  (IUNIT,  IER) 
where: 

ARGUMENT  DESCRIPTION  TYPE 


IUNIT  Logical  unit  number  of  photo  data  and  Integer 

directory  file 

IER  Returned  error  status  code  (see  Octal 

Appendix  4e  for  code  definitions) 
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ERROR  CODE  DEFINITIONS 


IER  (Octal)  =  error  status  code  of  the  I/O  function 


HR  TUHaTT - 

"  Explanation 

0 

Normal  I/O  completion 

1 

“End-of-file"  has  been  identified 

5 

Attempted  to  read  from  unassigned 
area  of  mass  storage 

20 

Write  attempted  on  read  only  file 

21 

Reference  made  to  unassigned  file 

22 

Write  attempted  beyond  assigned  area 
of  mass  storage 

25 

Maximum  length  of  input/output  array 
exceeded 

26 

FORTRAN  sentinels  are  unequal  -  read 
was  not  executed 

30 

Run  number  not  found  in  run 
di rectory 

31 

Variable  name  not  found  in  variable 
di rectory 

32 

Checksum  error 

33 

Run  number  not  last  entry  in  run 
di rectory 

34 

Update  in  place  aborted  due  to 
different  data  record  lengths 

35 

Invalid  FORTRAN  sentinel  -  word 
length  indicated  for  physical 
record  greater  than  249 

36 

Invalid  FORTRAN  sentinel  -  record 
number  Indicated  for  first 
physical  record  not  equal  to  1 
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